Extended droplet theory for aging in short-ranged spin glasses and a numerical 

examination 
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We analyze isothermal aging of a four dimensional Edwards-Anderson model in detail by Monte 
Carlo simulations. We analyze the d ata in the view of an extended version of the droplet theory 
proposed recently (oond-mat/020211C ) which is based on the original droplet theory plus conjectures 
on the anomalously soft droplets in the presence of domain walls. We found that the scaling laws 
including some fundamental predictions of the original droplet theory explain well our results. The 
results of our simulation strongly suggest the separation of the breaking of the time translational 
invariance and the fluctuation dissipation theorem in agreement with our scenario. 

PACS numbers: 



I 

C 

O 
o 



> 

m 
o 

(N 
O 

c3 



i 

a 

o 

o 



X 



I. INTRODUCTION 

Spin glasses exhibit characteristic slow dynamics below 
the spin-glass (SG) transition temperature T c . Recently 
aging phenomena n which have been known for a long time 
in glassy systemsEl, has attracted renew-ed|-wtejest both 
in experimental and theoretical studies&HotJ&ll A ma- 
jor theoretical progress was the development of the dy- 
namical mean-field theories (MFT) of spin glasses and 
related systemsEl. Most remarkably it found that the 
Field Cooled (FC) susceptibility xfc is larger than the 
equilibrium susceptibility xea which is related to the 
Edwards- Anderson (EA) order parameter gEA by the 
fluctuation dissipation theorem (FDT) as IcbTxea — 
1 — Sea- This finding corresponds to the well known 
experimental observationB that xfc is larger than the 
Zero Field Cooled (ZFC) susceptibility xzfc- In addi- 
tion, many new view points for the glassy dynamics were 
discovered subsequently such as the concept of effective 
temperaturetd'til. However, the MFT does not provide 
insights into what will become important in realistic fi- 
nite dimensional systems. Most seriously, thermally ac- 
tivated nucleation processes which are presumably im- 
portant in finite dimensional glassy systems cannot be 
captured at the mean-field level. 

Recently we proposed a refined scenariofc3 for the 
isotheijmaLraging based on the droplet theory for spin 
glassestjfiJO. We conjectured that the original idea of 
effective stiffness of droplets in the presence oi-frozen-in 
domain wall, introduced by Fisher and HuseE.3, can be 
extended to take into account anomalously soft droplet 
excitations which are as large as frozen-in extended de- 
fects, i. c. domain walls. This conjecture is partly moti- 
vated by the resplts of recent active studies of spin-glass 
models at T = Cm which revealed existence of anomalous 
low energy and large scale excitations. The anomalously 
soft droplets allow emergence of a new dynamical order 



parameter grj and the dynamical susceptibility xd asso- 
ciated with the former by FDT k^Txx> = 1 — <Zd- The 
dynamical order parameter go is expected to be smaller 
than the equilibrium EA order parameter ^ea which 
means that the dynamical susceptibility xd is larger 
than the equilibrium susceptibility xea- Consequently, 
our scenario implies a novel feature that breaking of the 
time translational invariance (TTI) and FDT separates 
asymptotically at large length (time) scales: the break- 
ing point of TTI converges to the equilibrium EA order 
parameter q^A while that of FDT converges to the new 
dynamical order parameter go smaller than qea- We will 
find that xfc — Xd so that our scenario also suggests 
Xfc > Xea- 

Both .-the original droplet theoryEj and our refined 
scenarioE3 predict scaling laws for the time dependent 
quantities measured in aging such as the magnetic auto- 
correlation function and the dynamical susceptibilities in 
terms of a time-dependent length scale L(t) which pre- 
sumably grows extremely slowly in a logarithmic fash- 
ion due to thermally activated processes. By now it is 
well understood that length scale that can be explored 
in practice is very much limited not only in numerical 
simulations (typically 1 — 10 lattice spacings) but also in 
real experiments (typically ~ 100 lattice spacings). The 
latter implies one must seriously take care possible pre- 
asymptotic behaviors to elucidate the desired asymptotic 
behavior associated with the putative T = glassy fixed 
point. To cope with such a complicated situation still 
in a controlled way, we examine the scaling theory by 
Monte Carlo simulations in two strokes. 

First, we examine the growth law of the dynamical 
length scale L[t) itself by directly measuring the spatial 
cohereHee-usittgi two real replicas. As realized in recent 
studiesOiislIj'E!!, the problem of crossover from critical 
to activated dynamics is the central issue here. Second, 
we examine the scaling properties of the time dependent 
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quantities of our interest by parameterizing the times us- 
ing the data of the time dependent length scale obtained 
by the separate simulation. The original droplet theory 
and our extended version provide some useful informa- 
tion of finite length correction terms to the asymptotic 
limit L — ► oo. The two-strokes (or parametric) strategy 
of the present paper, A^ichis already employed partially 
in the previous studiesBBBBEi allows us to cope with 
the mixture pre- asymptotic behaviors of different origins 
in a controlled way and far more advantages than usual 
approaches which try to examine the scaling laws in one 
stroke directly as a function of times blindly with many 
uncontrolled fitting parameters. 



In this paper we present a detailed study on the 
isothermal aging by Monte Carlo (MC) simulations on 
a four-dimensional (4d) EA Ising SG model. While the 
model system in 4d is somewhat non realistic, the advan- 
tage of studying the id Ising EA model is that impor- 
tant equilibrium properties concerning both the critical 
phenomena at T c and some essential scaling properties 
associated with the T — glassy fixed point, both of 
which will turn out to provide extremely useful infor- 
mation to study off-equilibrium dynamics, are far better 
known in Ad than in three dimensions (3d). In oder to 
take care of the critical fluctuations, we will use the well 
established information provided by the previous studies. 
Furthermore the value of stiffness exponent 8 is found to 
be considerably larger than that of 3d which implies eas- 
ier access to low temperature properties in Ad than in 
3d within limited length scales. Indeed, a recent anal- 
ysis of defect free-energy in Ad could clarify the antici- 
pated etpssover from critical regime to low temperature 
regimeEi In our analysis we employ the values of these 
parameters and fix them so that we are left with a few 
free parameters in our scaling analysis. 



The present paper is organized as follows. In the next 
section we introduce our model system studied. In sec- 
tion III, wc introduce the two time quantities used in 



Stion IV 



we ex- 



this paper and summarize some of their basic properties 
for the convenience of later sections. In 
plain our extended droplet scaling theory^ 



in a more self- 
contained and comprehensive manner recalling also t h e 
fundamental results of the original droplet theoryt^Edllla. 
In section we examine the growth law of the dynam- 
ical length L(t) by MC simulations. In section VI and 
VII, wc examine time dependent physical quantities by 
MC simulations and perform scaling analysis using the 
growth law L(t) obtained in section IV. A part of the 
resu lts w as already reported in Ref. |l7. Finally in sec- 
tion 
paper. 



II. MODEL AND SIMULATION METHOD 

We study the Ad Ising EA SG model, defined by the 
Hamiltonian 



H 



■ JijSiSj 



h^Si, 



(i) 



VIII, we present some discussions and conclude this 



where the sum runs over pairs of nearest neighbor sites. 
Ising variables are defined on a hypercubic lattice with 
periodic boundary conditions in all directions. The inter- 
actions are quenched random variables drawn with equal 
probability among ±J with J > 0. We will use J as the 
energy unit. The last term in the Hamiltonian represents 
the Zeeman energy with h being the strength of the ex- 
ternal uniform magnetic field. In this representation, the 
magnetic field h has the dimension of energy so that we 
will also measure it in the unit of J. In the simulations, 
we use k^T / J for the temperature scale and we set the 
Boltzmann constant — \ for simplicity. 

It has been well-established that a SG phase transition 
does occur at a finite temperature with strong ordering, 
namely a finite amplitude of SG order parametep-jn the 
ordered phase. Recent extensive MC studiesEffiS have 
estimated the critical temperature to be 2.0J. The crit- 
ical exponent v of the di verging coherence length is also 
obtained as v ~ 0.9 — LOnSEaO. Another important ex- 
ponent associated with T = glassy fixed point is the 
stiffness exponent 8 whose value is also obtained around 
0.7 by MC simulatiorc3 and ground state calculations. 
We note that the value of 8 ~ 0.7 in Ad is significantly 
larger than that of 3d Ising EA model, 8zd ~ 0.2a. This 
fact allows us to analyze the asymptotic behaviors rather 
easily than 3d case, which is one of our main reasons for 
investigating the Ad EA model. 

The simulation method is a standard single-spin-flip 
MC method using the two-sublattice dynamics with heat- 
bath transition probability. We define one Monte Carlo 
step (MCS) as N spin trials. We also use the multi-spin 
coding technique, which simulates 32 different systems 
independently at the same time on a 32-bit computer. 
The system sizes studied are L = 8, 16, 24 and 32 at T c . 
Below T C) we mainly study L = 24 systems, and L = 32 
in order to check finite-size effects. There is no significant 
difference between data of L = 24 and 32 at least within 
our time window (10 5 MCS). 



III. TWO-TIME QUANTITIES 

Experimentally, isothermal aging of spin glasses is in- 
vestigated by observing response of the system to an ap- 
plied external magnetic field. In the present paper, we 
study dynamical DC linear magnetic susceptibilities and 
their conjugate magnetic (spin) autocorrelation function 
during isothermal aging by Monte Carlo simulations. In 
the present section, we introduce the two time quantities 
and summarize some basic properties. 
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To mimic the experimental protocol of isothermal ag- 
ing, we consider that the configuration of the system is 
completely random at time t = and then start to re- 
lax in touch with a heat-bath at temperature T which 
is lower than the critical temperature T c . Thus cooling 
rate is infinitely fast. 

There are two standard protocols used in DC magne- 
tization measurements. In the so-called zero field cooling 
(ZFC) procedure, the system first evolves for a waiting 
time < w without an applied magnetic field then a small 
probing magnetic field of strength h is switched on. The 
growth of the induced magnetization is measured after- 
wards. In the measurement of the so-called thermorema- 
nent magnetization (TRM) , the system evolves under the 
applied magnetic field of strength h for the waiting time 
t w and then the field is cut-off. The decay of the mag- 
netization induced during the waiting time is measured 
afterwards. In our simulations, we measure the linear 
susceptibility at time t > t w as, 



X(Mw) 



i (Mm 

N h 



(2) 



where M(t) is the total magnetization measured at time 
t with N being the number of spins. The total magneti- 
zation is given by sum 



M(i) 



(3) 



where Si(t) is the Ising spin variable at site i at time t 
after the quench. Correspondingly, we measure the mag- 
netic (spin) autocorrelation function, 



C(Mw) - ^(M(t)M(t w )) = £>«(W*")> 



(4) 



The last equation holds for our model with no ferromag- 
netic or anti-ferromagnetic bias so that the sum of the 
cross terms of i ^= j vanishes as 1/ \/~N — > in the ther- 
modynamic limit N — ► oo. Furthermore in our Ising 
model, the spins are normalized Sf — 1 which yields, 



c(*,*) = i£<Si(i) 2 ) = i 



(5) 



for any t. 

In the above equations, bracket (. . .) means to take an 
average over different realizations of initial conditions, 
thermal noises and realization of random exchange cou- 
plings of the system. However, the above quantities are 
presumably self-averaging for thermodynamically large 
systems N — > oo. 

If linear response holds, the dynamical linear suscep- 
tibilities measured in the ZFC and TRM procedure can 
be written as, 



dt'R(t,t') 



(7) 



where R(t, t') is the magnetic linear response function. 
The latter is define as, 



v ' ; TV Sh^o Sh(t') 



t > t' 



(8) 



where S(M(t)} is the induced magnetization at time t by 
an infinitesimal probing pulse field 5h(t') applied only at 
time t'(< t). From these, it follows that the sum of the 
two susceptibilities 

Xzfc(Mw) + Xtrm(Mw) = / dt'R(t,t')= X ZFc(t,0) 



(9) 

becomes independent of the waiting time i w and only a 
functioH|Of the total time t elapsed after the temperature- 
quenchQ'tl One can use this sum .ride as a criterion to 
check linearity of measurementsoEjtEl. 

Let us briefly discuss the implication of the sum rule 
(^) combined with the following very mild assumptions. 
First, the ZFC linear susceptibility lim tw ^o Xzfc(£, ^w) 
increases with t but saturates since it is bounded from 
above. Let us define in particular the limit, 



Xfc = Am Xzfc(*,0). 



(10) 



We call the latter as Field Cooled (FC) susceptibility 
because it is more or less similar to what is called FC 
susceptibility as we discuss below. _S|econd, it is natural 
to assume weak long term memoryjtj, 



lim xtrm(Mw) 

t — »OG 







(11) 



for any large but finite t w . The latter simply means that 
TRM should relax down to zero for any large but finite 
waiting time t w during which the magnetic field is ap- 
plied. Then we find using above assumptions in the sum 
rule (0), 



lim xzfc(Mw) = lim Xzfc(*, 0) 

t — >oo t—*oo 



Xfc 



(12) 



for any large but finite t w . 

Let us explain why we call ( |l0| ) as the FC susceptibility. 
Usually the FC magnetization is measured by cooling 
down the temperature with a certain cooling rate from 
above T c down to a target temperature T below T c with 
the magnetic field h being applied. Suppose that it takes 
time e to cool down the temperature (typically of order 
30 sec) and the target temperature is reached at time 
t = 0. Then the linear susceptibility measured in this 
protocol can be expressed as, 

,o 

M FC (t)/h= / dt'R(t,t')+ X ZFc(t,0). (13) 



Xtrm(£> t w ) 



dt'R(t,t') 



(6) 



The above expression is formally valid as far as lin- 
ear response holds. Note that the response function 
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in the first term is defined with respect to the partic- 
ular schedule of the temperature changes. The contri- 
bution of the first term decreases with time t because 
of the weak long term memory property (|fl|). Thus 
in the limit t — > oo the susceptibility Mpc(t)/h con- 
verges to the FC susceptibility of our definition in (10), 
lim^oo M FC (t)/h = limt^oo xzFc(t,0) = xfc- The ap- 
proach to the limit may well be slow. Experimental ob- 
servations (See for instance Fig 13. of ||) show that the 
correction term to the asymptotic limit relaxes slowly 
but the amplitude is very small such that it can be made 
much less than 1% of the asymptotic value well within 
the experimental time window. 

In the present paper, we do not discuss the possible 
effects of finite cooling rates and assume the idealized 
temperature quench e = 0. In this idealized situation, 
the first term in ( |l3| ) is absent and the TRM becomes 
equivalent to the so-called isothermal remanent magne- 
tization (IRM). 

Another interesting limit is to consider t w — > oo first 
with fixed time separation r = t — t w . In this limit, one 
expects to find equilibrium (stationary) response, 

(14) 



Xeq(l~) = lim XZFCO +*w,*w)i 



which only depends on the time separation r. The equi- 
librium susceptibility xea is defined as, 



XEA = lim Xeq(T). 



(15) 



The last static susceptibility xea is more or less close to 
what is called the ZFC magnetization (divided by h). 

A very important issue is then whether the two suscep- 
tibilities xea and Xfc the same or different. As we 
noted in the introduction this is intimately related with 
the fundamental experimental observation in spin-glass 
systemsEI, namely xfc > Xzfc- One of the most remark- 
able finding of the dynamical mean-field theorytl is that 
indeed an inequality 



Xfc > Xea 



(16) 



holds with xfc defined in ([H]) and xea defined in 
(|l5|). The difference is due to anomalous contribution 
of the slowly relaxing, aging part of the response func- 
tion i?Ii 1( £A_.On the other hand, the conventional droplet 
theoryE3oll3 was understood!] to predict xfc = Xea , 
i. e. no anomaly. As we explain later in section IV D, 
our extended droplet theorycHI predicts (16) with xfc 



being identified with the dynamical susceptibility xd 
associated with the noble dynamical order parameter 

qu{< <7ea) as k B TxF> = 1- qu- 

Another important issue is to what extent FDT holds 
in aging systems. In our present context FDT reads as, 

R(t,t') = ^d t ,C(t,t') (17) 

where C(t, t') is the autocorrelation function. It becomes 
after integration over time J t dt' ■ ■ • , 



Since this must hold precisely in equilibrium, the equi- 
librium limit of ZFC linear susceptibility ( |l4| ) must be 
related to that of the spin autocorrelation function as, 



Txc q (T) = 1 - C eq (r) 



where 



a 



cq 



(r) = lim C(t + t Wl t w ), 



(19) 



(20) 



is the spin autocorrelation function in the equilibrium. 
In the static limit t — > oo ( |T9| ) becomes the static FDT, 

^bTxea = 1 - <?ea, (21) 
where qea is the static EA order parameter defined as, 



<7ea 



lim Ceq(r). 



(22) 



Except for the ideal equilibrium limit, FDT (|l^) is not 
guaranteed in general. However, it was realized recently 
by Dean, Cugliandolo and Kurchan that possible ampli- 
tude of the violation of the FDT 0, 



I(t, t w ) = l- C(t, t w ) - k B TxzFc(t, t w ) 



(23) 



is bounded from above by the entropy production rateEl 
The bound implies even for very slowly relaxing sys- 
tems in which entropy production rate becomes small, 
the FDT (|l8|) should hold between spontaneous thermal 
fluctuations and linear responses at least for short enough 
time scales r. 

Let us consider asymptotic limit t w — > oo of the two 
time quantities with fixed value of the autocorrelation 
function C = C(t + t w ,t w ). First let us note that 
linv^oo C cq (T) = <7ea implies the time translational in- 
variancc (TTI) is strongly broken at < C < qEA, i- c. 
the two time quantities cannot be a function of only 
the time separation in this regime. Concerning the in- 
tegral FDT violation (p3|), it becomes a function of C, 
i. e. limt w _>co I(t, t w ) = 1(C). In conventional cases 
I (Sea < C < 1) = while 1(0 < C < 9e a) > bjua 
a non-trivial function of C by the dynamical MFTtrEj 
and 7(CL^i-|C. < <7ea) = ?ea — C in usual coarsening 
systemscilE§E3. On the other hand, our scenaricO sug- 
gests I(q D < C < 1) = and 1(0 < C < q D ) = qr> - C. 
Here q^> is the new dynamical order parameter which is 
smaller than qea- Thus the breaking points of TTI and 
FDT take place separately at qea and qn respectively 
in our scenario while both of them take place simultane- 
ously at. gEAiiirf the- dynamical MFT and usual coarsening 
systemsoBEj. 



IV. THEORETICAL BACKGROUND 

In this section we discuss our extended droplet theory 
sketched in Ref. 12 concerning isothermal aging which 



FDT 



1 - C(t,t w ) = k B TxzFc(t,t w ). 



(18) 



will be our basis to analyze the data of Monte Carlo simu- 
lations in later sections. As we noted in the introduction, 
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we pay a special attention to the idea of so-called effec- 
tive stiffness of droplet excitations in the presence of do- 
main walls which are present as extended defects during 
isothermal aging. For clarity, we will try to present this 
section in a self-contained fashion including , summaries 
of the results of the original droplet theoryoEj which 
are almost fully included in our scenario. To simplify 
notations, we consider systems of N Ising spins Si = ±1 
(i = 1, . . . , N) in a d-dimensional space coupled by short- 
ranged interactions of energy scale J with random signs 
with no ferromagnetic or anti- ferromagnetic bias. 

A. Basics 

Let vts, recall briefly the starting point of the droplet 
theory.E-3 It assumes that thermodynamic states of (Ising) 
spin-glass phases consist of a pair of pure states of an in- 
finite system which are related by global spin inversion. 
At low but finite temperature in the spin-glass phase, an 
equilibrium state can be considered as made of a ground 
state, say T or its global spin inversion T, plus thermally 
activated droplet excitations of various sizes taking place 
on top of r. In simple systems such as ferromagneis 
droplet excitations exist but play a rather limited -role ia. 
An essential finding of the original droplet theoryldo is 
that the temperature is dangerously irrelevant in spin- 
glass phases because of strong impacts of thermally acti- 
vated droplets. 

A droplet at a given length scale L is supposed to be a 
compact cluster of spins with a volume L d with d being 
the dimension of the space and a surface volume L f with 
df being the (fractal) surface dimension. The typical 
value of its excitation gap F^ p with respect to the ground 
state is supposed to scales as 

p ~ T(L/L ) e , (24) 

where the exponent 8 > is the stiffness exponent, T 
is the stiffness constant and Lq is a microscopic length 
scale. The excitation gap is however broadly distributed 
with the typical value given above. The probability dis- 
tribution of the free-energy gap is expected to follow a 
universal scaling form, 

P (F L )dF L = p(F L /F^)dF L /F^ (25) 

with non-vanishing amplitude at the origin 

p(0) > (26) 

which allows marginal droplets. This property p(0) > 
allows dangerously active droplets which will respond 
to arbitrarily weak perturbations so that they play ex- 
tremely important roles as the Goldstone modes: they 
dominate spontaneous thermal fluctuations and linear re- 
sponses. 

Dynamically, the excitation of such a cluster of spins 
is supposed to happen only by thermal activated process. 



The typical value of free-energy barrier B^ p to flip the 
cluster of spins is supposed to scale as 

~ A(L/Lo)* (27) 

where A is a characteristic free-energy scale of the barri- 
ers. The Arrhenius law implies that a droplet of length 
scale L(t), 

£(*)~(4-k(*Ai>)) . (28) 

can be activated within a time scale of t. Here To is a 
certain unit time scale for the activated processes. Let 
us call this time dependent length scale as the dynamical 
length scale. 

B. Crossover from Critical to Low Temperature 
Regime 

In practice, it is necessary to take into account of crit- 
ical fluctuations near T c . Even at T < T c , the length 
scales shorter than the coherence length of the critical 
fluctuation 

£_~L |l-r/T c |-" (29) 

should be dominated by critical fluctuations. Asymptotic 
low temperature properties should appear only at larger 
length scales. 

Correspondingly, we expect two typical stages in the 
dynamical length scale L(t) as shown in Fig. Q. One is a 
critical dynamics associated with critical slowing down in 
time range r (T) 3> t >• t , where L(t) follows a power 
law with the dynamical critical exponent z, 

L(t) = / (iAo) 1/z . (30) 

Here Iq is a microscopic length scale which is of order 
1 lattice distance in EA models and to is a microscopic 
time scale which is typically to ~ 10~ 12 — 1CP 13 (sec) in 
real spin systems while it is 1 Monte Carlo Step (MCS) 
in usual heat-bath Monte Carlo simulations. This for- 
mula means that even at T < T c , L(t) behaves like the 
critical power-law at short time-length scales. On the 
other hand, the intrinsic low-temperature dynamics as- 
sociated with the T — glassy fixed point ( p8| ) should 
occur only at larger time scales beyond a certain crossover 
time t > t (T). 

The crossover length Lo(T) would be determined by 
a comparison between free-energy barrier and thermal 
energy, 

k B T~A(T)(L (T)/l )l>, (31) 

where the characteristic free-energy scale A(T) behaves 
like A(T) = J|l - T/T c \^ near TjB- It leads to the 
temperature dependence of Lq(T), 

L (T) = l (T/J) 1 ^\l-T/T c \-\ (32) 
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Time or 

spatial length scale 




r < 77/ 

FIG. 1: Schematic picture of length scale. 

The corresponding 



which is essentially equivalent to (|2 
crossover time tq(T) is given by 

r (T) = t (L (T)/l ) 



(33) 



We obtain the singular part of the crossover time at T, 
To(T) ~ |1 — T/T c \~ zu , as expected from a critical scaling 
theory. Consequently, the scaling formula of the growth 
law which describes the whole crossover from the critical 
dynamics at t <C tq(T) and the activated dynamics at 
t > t (T) is given by 



L(t)/L (T) = L(t/T Q (T)), 



where 



L(x) 



x 1/z (x<Sl), 
log^ix) (z»l). 



(34) 



(35) 



One should note that the crossover could be very gradual 
and functional form of the intermediate regime (which 
will dominate realistic time ranges in simulations and 
experiments) can have very complicated expression which 
is not obvious. This crossover in the growth law of the 
dynamical length scale is numerically examined in section 
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The importance of the crossover from critical to low 
tepaperature behavior has been pointed out by Bokil et 
alEZI concerning some static properties of low-temperature 
SG phase. Let us note that the above analysis of the 
crossover from the critical to activated dynamics is a di- 
rect (dynamical) analogaie of the analysis of defect free- 
energy in Ad EA modelo by one of us (KH) . In the lat- 
ter study, the defect free-energy was found to become 
a universal constant in the limit L <C £- and grows as 
T{T){L/L ) e and T(T) - J|l - T/T c \ e » with 6 = 0.82, 
v = 0.93 at L > 



C. Domain Growth 

During isothermal aging up to a time t after quench, 
domains with the mean size L(t) separating different pure 



states grow up by coarsening domain walls of smaller 
length scales. The droplet theory proposed scaling prop- 
erties of time-dependent physical quantities in terms of 
time-dependent mean domain size L(t). We discuss the 
scaling by L(t) of one-time quantities with time elapse 
after the quench, which examine later by MC simula- 
tions. One example is time development of energy per 
spin defined by 



(36) 



Relaxation of the energy per spin is expected to be due 
to relaxation of excessive energy associated with domain 
walls. Thus it is expected to decrease asE3, 



e(t) 



--eq 



In 



(37) 



where T' is a temperature dependent parameter. 

Another interesting one-time quantity is domain-wall 
density p s (t) in which the morphology of the domain 
with the fractal surface dimension df > d — 1 appears. 
In coarsening dynamics the density decreases with time 
during isothermal aging. In simple systems such as a 
ferromagnet where the domain wall becomes flat at suf- 
ficiently low temperature, the density of domain wall is 
proportional to the inverse of the mean size 1/L(t). How- 
ever, it could be rough with the fractal dimension in spin 
glasses because of the disorder and frustration so that we 
expect 



Ps(t) 



lo 



(38) 



In MC simulations, two replicas with identical interaction 
bonds are updated independently and the domain-wall 
density is calculated afl 

Psit) = \ (l -±J2*l a \t)*f(t)*^(t)a?\t)) , 

\ (ij) / 

(39) 

where Nb is the number of bonds and the suffixes a and 
(3 denote replica indices. Here, following Ref. ^l] we take 
short time average of each spin as 



(t) = sgn 




(40) 



This procedure in ( f4fj| ) means that smaller fluctuations 
associated with small droplets are eliminated and the 
coarsening domain walls arc emphasized. Without taking 
such average over time, the density is similar to so-called 
link overlap function, which does not vanish in the long 
time limit. 
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D. Domain Walls and Soft Droplets 



Two Length Quantities 



Following Ref. 15, the whole process of isothermal ag- 



ing may be divided into epochs such that the typical sep- 
aration between domain walls (hereafter simply denoted 
as domain size) Lo, aLg, a 2 Lo, . . . , a n L>o, . . . where a > 1. 
At each epoch, droplets of various sizes up to that of the 
domain size can be thermally activated or polarized by 
the magnetic field. The two time quantities we intro- 
duced in section [II namely the autocorrelation function 
C(t + t w ,t w ) defined in (M) and the linear susceptibility 
x(t + t w ,t w ) defined in (g) probe such thermal fluctu- 
ations and linear responses of droplets smaller than the 
size of the domain if the time separation is limited such 
that L(t) <L(t w ). 

In the previous workli3, we extended the idea of Fisher 
and Husdla who noticed that droplet excitations can be 
softened in the presence of a frozen-in defect or domain 
wall compared with ideal equilibrium with no extended 
defects. This is because droplets which touch the defects 
can reduce the excitation gap compared with those in 
equilibrium. This effect will have very important impacts 
on the two time quantities. 

Let us consider a system with a frozen-in defect of size 
R: a large droplet of size R is flipped with respect to T, 
which is a ground state of an infinite system, and then it 
is frozen. The typical free-energy gap F^j, of a smaller 
droplet of size L in the interior of the frozen-in defect is 
expected to scale as 



FZ P R = T eS [L/R](L/L ) 



L < R 



(41) 



To explain the consequence of our conjecture above in- 
troduced, let us consider thermal fluctuations and mag- 
netic linear responses of droplet excitations of size L 
enclosed in a compact region made by a frozen-in de- 
fect at scale R which has iV<j oc (R/Lo) d spins— To 
this end, let us construct a toy droplet modelEJ de- 
fined on logarithmically separated shells of length scales 
Lfc/Lo = ak < R/Lq with a > 1 and < k < n-x < n\ 



where R/L = a ni and L m /L 



For each shell 



an optimal droplet is assigned whose free-energy gap is 
minimized within the shell. Then F L y j, will be of order 
Fq at the shell k = m so that we have 



ptyp 

L.R 



ptyp 

L.R 



(1 



F S k , 



(45) 



For simplicity, droplets at different scales are assumed to 
be independent from each other. 

At the shell L/Lq = a k , the system is decomposed into 
compact cells of volume L d such that each cells represents 
a droplet of size L. The number density of droplets per 
spin associated with the shell scales as, 



In a. 



(46) 



Each droplet excitation will induce a random change of 
the magnetization of order 



M L ~ m^J{L/L ) d , 



(47) 



with T c g[L/R] being an effective stiffness which is only 
a function of the ratio y — L/R. _ 
For y <C 1, T G ff [y] will decrease with y aaij 



T e ff[y]/T = 1 - c v y 



d-8 



for 



(42) 



Here T is the original stiffness constant T = T e ff(0). A 
basic conjecture, on which our new scenario based, is that 
at the other limit y ~ 1 the effective stiffness vanishes as, 



TeffM/T ~ (i - y y 



(43) 



with < a < 1 being an unknown exponent, and that 
the lower bound for F^^ should be of order J, say F . 

We assume that the probability distribution of the free- 
energy gap Fl,r is broad and obey the same functional 
form as ( p5|) where F l L yv should be replaced by Ff^L. 
Then the scaling form of the distribution of the gap (Eq) 
and non-zero amplitude of gap-less droplets p(0) > ( pq ) 
implies that the probability that a gap is smaller than a 
certain threshold 5U(<g. F L yp ) scales as, 



where m is the average magnetic moment within a vol- 
ume of L . 

The thermal fluctuation can be measured by an order 
parameter 



<1 



2 

%i i 



(48) 



where the sum runs over sites in the interior of the frozen- 
in defect. Here we have put the over line 7TT which means 
to take average over different realization of the random- 
ness. In the absence of any droplet excitations, q = 1 
holds due to the normalization of the Ising spins. A spon- 
taneous thermal fluctuation of a droplet will take place 
if its excitation gap happens to be smaller than the ther- 
mal energy k^T. The probability of the latter is found 
to be proportional to p(0)fc B T / 'F^ due to @. So the 
reduction from 1 due to a droplet excitation at scale L is 
of order Mlp(0)(k B T '/ 'F^) . Then the reduction of the 
spin autocorrelation by droplet excitations at scale L is 
estimated as, 



Prob(F L , i j(<5[/) 



(44) 



Since the probability is linear mthj-^U, these active 
droplet are dangerously irrelevantl 13 rl 1 'ii 15 l: arbitrary small 
perturbation 5U may trigger a droplet excitation. 



5g L ~Mlm^(j£) In a (49) 

where the last factor is due to the number of droplets per 
spin given in ([46|). 
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The magnetic linear response by weak external mag- 
netic field h is measured by a linear susceptibility 



X =hvaN-^{Si)h/h, 



(50) 



where ^2i(Si)h is the induced magnetization by the held. 
Since a droplet excitation can induce a magnetization of 
order Mr, given in ( [47| ) with random signs, it can gain 
a Zeeman energy of order SU l ~ Mr,5h by responding 
to the held. The probability that the droplet excitation 
takes place is then proportional to the probability that 
the gain by the Zeeman energy is larger than the exci- 
tation gap of the droplet which can be found using (44). 
Then the expectation value of the magnetic moment in- 
duced by droplet excitations at scale L is estimated as, 



hS XL ~ M L p(0) 



hM L 
FT 



In a. 



(51) 



Now summing over contributions at different length 
scales < k < ri2, we obtain the total reduction of the 
order parameter from 1 and the linear susceptibility as 
using ©, © and @ by, 

1 - q(L m , R) = fc B Tx(i m , R) 



p(0)m 



2 ^bT 

fc=o F l,rQ ~ ^"i ) + F oSk,m 



p{Q)m 2 k B T x 



Lm dL 
T 



L 



l-A (ln(L/i?)) , A Q (ln(L/i?)) 



ptyp 
L,R 



Fn 



(52) 



In the last equation, the sum X)fc=o * s replaced by an 
integral J^" 1 dL/L and A a (z) is a pseudo ^-function of 
width lncM Note that FDT is satisfied between @ 
and @. 

In the following let us call the two length quantities 
q{L, R) and x(L, R) as the generalized order parameter 
and the generalized linear susceptibility, respectively. We 
will associate these two-length quantities with the two- 
tame quantities measured in aging experiments. 



2. Edwards-Anderson Order Parameter 

For clarity, let us consider the equilibrium limit where 
there is no extended defects which can be realized by 
taking R — > oo hrst. In the latter iipiit we recover the 
result of the original droplet theorytil, 

p(0)m 2 k B T , , 

^hm q(L, R) = q EA + c , (53) 



with 



/oo 
dyy' 1 ' 



(54) 



and the Edwards-Anderson (EA) order parameter de- 
fined in (E2|) evaluated as, 



(55) 



k B T 

q E A = hm lim q(L,R) = 1 - cp(0)m 2 — — . 

L — >oo R — >oo 1 



The associated equilibrium susceptibility xea is dehned 
as ( pl| ) fcB^XEA = 1 — qEA- The above expressions be- 
come useful when we consider equilibrium dynamics. 



3. Dynamical Order Parameter 

It is useful to consider asymptotic behavior at large 
sizes R/L ^> 1 with the ratio x — L m /R being hxed. 
We obtain for < x < 1, 

< f? m i P(0)™ 2 kBT p{Q)m 2 k B T ^ 
q(xR,R) = q E A+ ^ /U/T w A{x) 6a(ma;) 



T(i?/L ) 6 



with O a (z) being a pseudo step-function of width hn 
and 



A{x) 



dyy 1_ 



dyy- l - e {T/T cS [y]{\ - A a (lny)) - 1). (57) 



Note that the second integral coituerges because of ( |42| ) 
and the inequality 9 < (d — 1)/2l3. Thus as far as < 
x < 1, the order parameter converges to the equilibrium 
EA order parameter (jea evaluated in (55|). 

In the intriguing case x ~ 1, A{x) will remain hnite as 
far as < a < 1. At x ~ 1 the last term of (p6[), which 
is due to the anomalously soft droplets, contributes and 
we obtain, 



q(R, R) = q D + A(l) 



p{0)m 2 k B T 



T(R/L )o ' 

where we have dehned the dynamical order parameter 



(58) 



qo 



lim q(R,R) 

R — >oo 



qEA - p(0)m 



k B T 



2 K B 



F 



(59) 



Naturally, we can dehne the associated dynamical linear 
susceptibility xd as 



k B TxD = l-qr>- 
The above results imply 

qu < qEA 

and 

XD > XEA 



(60) 



(61) 



(62) 



As we discuss below go and xd play significantly im- 
portant roles in the dynamical observables of aging. We 
will see that the field cooled susceptibility xfc defined in 
( |l0| ) is equal to the dynamical susceptibility xd dehned 
above. Thus the inequality (|6^) implies the anticipated 
inequality xfc > Xea given in (pi). 
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4- Two-time and Two-length Quantities 



1. Equilibrium Limit 



We can now associate the two length quantities dis- 
cussed above with the two time quantities for short time 
separations L(t) < L(t w ). In the latter regime, the au- 
tocorrelation function C(r + t w ,t w ) and the ZFC linear 
susceptibility Xzic{ T + t w ,t w ) measure respectively ther- 
mal fluctuations and linear responses of droplets smaller 
than the size of the domain L(t w ). 

The autocorrelation function at a given time separa- 
tion t probes thermal fluctuations of droplets as large as 
L(t) in the presence of domain walls of size R(t w ) so that 
we expect 



C(T + t w ,t w ) = q(L{T),L(t w )) 1 



(63) 



where q(L, L') is the generalized order parameter given 
in (52). Similarly, the ZFC susceptibility probes the po- 
larization of droplets as large as L{t) in the presence of 
domain walls of size R(t w ) we expect, 



kBTxzFc(i~ +t w ,t w ) 



k B T X (L( T ),L(t w )) 
l-C(r + t w ,t w ), (64) 



where x(^j L') is the generalized susceptibility given in 

Let us note that the generalized order parameter 
q(L, L') and the generalized susceptibility x(£, V) are 
defined via disorder averaging of many different realiza- 
tion of small systems of size L'. Here it is important to 
recall that at any finite time, a macroscopic system will 
contain macroscopic number of domains no matter how 
large their size L' = L(t w ) is. Thus we can safely evaluate 
these two time quantities, which are macroscopic quan- 
tities, by the disorder-averaged two-length quantities. 

Let us emphasize that the FDT (|l^) is satisfied while 
TTI (Time Translational Invariance) is broken in the 
above two length/time quantities for L(t) < L(t w ). The 
latter is due to the fact that the autocorrelation function 
and the susceptibility depends on two length/time, i.e. 
not only L(t) but also on L(t w ). 

For larger time separations L(t) > L(t w ), we take into 
account decay of memory due to coarsening of domain 
walls following the original droplet theorylij. 



E. Scaling Properties of Two-Time Quantities 

Using the results of the previous section, we now dis- 
cuss scaling properties of the two time quantities at dif- 
ferent regimes in detail. Basically we expect three dis- 
tinct regimes i) quasi-equilibrium regime L(t) < L(t w ), 
ii) crossover regime L(t) ~ L(t w ) and iii) aging regime 
L(r) > L(t w ). In the following we first consider the ideal 
equilibrium limit for clarity and then the three regimes 
subsequently. 



Let us consider first the equilibrium limits of the auto- 
correlation function and the ZFC susceptibility which are 
obtained by taking the limit L( t w ) — > oo with fixed L(t) 
in (|ipt) , (|63|) and 
findsE 



From ( [V D 2 ) one immediately 



C cq (r) = q EA + cp(0)m 2 (k B T/T)(L(T)/L y 



(65) 



and 



Xcq(r) = lim xzfc(t + t w ,t w ) = XEA - 



p(0)m 2 



T(L(r)/L ) e 
(66) 

In (|6^) q-EA is the EA order parameter defined in ( j22] ) 
and evaluated in (|55"|). The numerical constant c is given 
in (|54"|). Correspondingly xea is the equilibrium suscep- 
tibility related to (Jea as k B TxEA = 1 — 9ea as in (^l|). 



2. Quasi-equilibrium regime 

Now we discuss the quasi-equilibrium regime L(t) < 
L(t w ). The autocorrelation function and ZFC linear sus- 
ceptibility can be evaluated using 



063) and fl6j). 

For simplicity, here we assume to stay deep in the quasi- 
equilibrium regime L(j) -C L(i w ). Then using the scaling 
property of the effective stiffness (ff2|), one obtainscS, 



c'p(0)m 2 



t -4" t w , iw) = 

knT 



a 



eq(r) 

L(r) 



T(L(t)/L ) 6 \L(U,) 



d-i 



, (67) 



with C eq (r) being the equilibrium part given in ( |65| ) and 

d = c v dyy 2 ^ d ~ 1 / 2 ^ \ The 2nd term is the weak non- 
equilibrium correction term due to the weak softening of 
small droplets which gives rise to some weak waiting time 
dependences. 

The above formula combined with the formula for the 
equilibrium part (|6^) yields a systematic extrapolation 
scheme to determine the EA order parameter ^ea- Such 
an extrapolation was al ready demonstrated in our previ- 
ous papertll. In section VII C , we perform such an anal- 
ysis numerically. 

Similarly, the ZFC susceptibility xzfc( t + t w , t w ) is 
evaluated as, 



kBTxzFc(t 
+c'p(0)m 2 



1 k^Tx 
L{r) 

T(L(T)/L ) e \L{t w ) 



t -\- t w , £ w ) 
k B T 



(r) 



eq 
d-8 



(68) 



|) and 



with Xcq( T ) being the equilibrium part given in 
Xea is the equilibrium susceptibility. 

The quasi-equilibrium regime is most relevant for the 
measurements of the relaxation of AC susceptibilities 
during isothermal aging. Recently a related scaling anal- 
ysis was performed in an experimentlij. Let us note 
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that previous experimental analysis of spontaneous ther- 
mal fluctuation of the magnetization and AC susceptibil- 
ity have already confirmed that FDT holds for quasi- 
cquilibj^ um , regime while non-stationarity is observed 
clearlyc 



3. Crossover Regime 

Next let us consider the crossover regime L(t) ~ L(t w ). 
Here we have to consider also the contribution of anoma- 
lously soft droplets of sizes as large as that of the domain. 

From ( J5S| ) the spin autocorrelation function is obtained 
immediately as, 



C(r + i w , fw)|z,( T )/L(t w )~i 



q D +p(0)m 2 A{l) 



T(L(t w )/L ) s 



(69) 



where is the dynamical order parameter defined in 
(|59|). Similarly the ZFC susceptibility is obtained as, 



&bTXZFc(£w + t w , t w )|i(r)/L(t w )~l 



k B T X D - p{Q)m z A{\) 



T(L(t w )/L o y 



(70) 



where \t> is the dynamical susceptibility defined in (|6C 
As we discuss in section IV E 6, the change from the 



quasi-equilibrium regime and the crossover regime be- 
comes very abrupt as function of x = L(t)/ L(t w ) in the 
asymptotic limit L(t w ) — > oo (See Fig. ||). This novel 
feature deserves to be studied in simulations and experi- 
ments. A convenient measure is the modified relaxation 
rate function S mo< i(x,t w ) we introduce in (|9l|). 



4- Aging regime : autocorrelation function 

Let us now consider aging regime L(t) ~ L(t) > L(i w ) 
where we need to consider growth of the domains ex- 
plicitly. For simplicity, we used the notion of epochs 
mentioned previouslylij. The n-th epoch spans logarith- 
mically separated time scales between £ n _i and t n such 



that L(t„_i) 



Lq and L{t n ) = a n Lo with a > 1. 



At each epoch some of the smaller domains are elimi- 
nated so that domain walls are coarsened. Thus a given 
site may or may not belong to the same domain (r or 
f) at two different epochs. The probability P s (i?i,i?2) 
that a given site belongs to the same domain at the two 
different epochs characterized by the sizes of the domains 
i?i and i?2 presumably becomes a function of the ratio 
of the size of the domains, i. e. 



(71) 



reflectin 
processc 



-the self-similar nature of the domain growth 
!i This function describes the slow decay of 



memory by the domain growth and plays the central role 
in the following. 

For simplicity, let us neglect thermally active droplet 
excitations in the interior of the domains and consider 
only the domain growth itself for the moment. Then the 
auto-correlation function between two epochs such that 
size of domains are R\ = L(t) and i?2 = L(t w ) becomes 



C(t = T + t w ,t w ) = C 



L(t) 
L(t w ) 



with 



C(x) = 2P s (x) - 1. 



(72) 



(73) 



In the limit x — > 1 , the scaling function should converge 
as 



lim C(x) = 1, 

SC— 3-1+ 



(74) 



because the normalization of the Ising spins. In the other 
limit, the scaling function is expected to behave asymp- 
totically as 

C{x) - x~ x (x > 1) (75) 

with A being a non-equilibrium exponent0El 

d/2 < A < d. (76) 

By taking into account thermally active droplet exci- 
tations, we expect the scaling form of the spin autocor- 
relation function as 



C(t,t w )~q u C 



L{t) 
L(t w ) 



(77) 



where qd is the dynamical order parameter . introd uced 
in ( |59| ) . Note that in usual coarsening processE§E3ci, the 
amplitude is given by EA order parameter g&A which was 
also assumed in the original droplet theorycJ. However, 
we expect the dynamical order parameter qD is more nat- 
ural because of the anomalously soft droplets which are 
as large as domains. The above scaling form ( f77|) should 
be manifested in the asymptotic limit L(t w ) — > oo with 
the ratio x — L(t) / L(t w ) fixed to a certain value larger 
than 1. Note that the normalization (|74|) allows match- 
ing with the crossover regime discussed in the previous 
section in the asymptotic limit L(t w ) — + oo. 

At finite time scales, we expect some correction terms 
because of the following reasons. First let us recall that 
the asymptotic amplitude of the oder parameter is at- 
tained only by integrating out contributions of droplet 
excitations up to infinitely large ones. At finite length 
scales, the amplitude of the order parameter should be 
larger than go due to the algebraic correction term such 
as the second term in (p8|) or (69). Then the factor qn in 
the scaling form ( |77| ) should be replaced by some time- 
dependent factor at finite time scales. Second, we also 
expect some additive correction terms since with some 
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probability there will be some region where domain walls 
don't pass through so that the dynamics due to droplets 
in the interior of the domains continues. Such quasi- 
cquilibrium corrections will be additive as considered in 
dynamical MFTu. However, we don't know how to col- 
laborate both the multiplicative and additive correction 
terms simultaneously and we will leave the problem of 
correction terms for the autocorrelation function in the 
aging regime for future studies. 



5. Aging regime : linear susceptibilities 

Finally let us consider the linear response to magnetic 
field during the domain growth (aging). Suppose that 
magnetic field h is applied only during a certain epoch 
where the size of the domain is L n / L Q — a n with n be- 
ing a certain positive integer. Then let us consider the 
magnetization SM n (t) measured at some time t > t n —\ 
during and after the epoch, 



SM n (t) 



-i(t,t„) 



dt'R(t, t') 



(78) 



t n -i 



where R(t, t') is the response function defined in (Q). 

During the epoch, droplets over the length scales from 
L up to L(t = t — t n _i) can be polarized by the field. 
Importantly the size of the domain can be considered 
as frozen in time to the value L n during this epoch so 
that the linear response is the same as that in the quasi- 
equihbrium/ crossover regimes. Thus the magnetization 
measured within the same epoch will be 

5M n (t) ~ hx(L(r),L n ) 



hp(0) 
for 



W dL 



rrr 



L F typ 

L,L„ 



L 

t n -! <t<t 



(79) 



where we used the generalized susceptibility given in (52). 
This part due to the droplets satisfies the FDT (|is|). 

When the field is cut off at time t n , de-polarization of 
the droplets will start. Let us now follow the argument 
of Fisher and HuseE£l, and first consider a certain early 
stage of the n + 1-th epoch, say up to time i+ slightly 
after t n , such that further domain growth still does not 
proceed appreciably. Note that the switching-off of the 
field is equivalent to adding additional field of the oppo- 
site sign —h. The latter will induce additional negative 
magnetization whose amplitude grow like ([79|) with r be- 
ing understood as the time after the field change. Up to 
the time i+ most of the magnetizations will be canceled. 
However some residual magnetization of order, 



6M n (t+) ~ c"p(0) 



hm 2 



r(L n /L y 



(80) 



may be left behind. Here c" is a certain numerical con- 
stant. As domain growth proceeds further, the residual 



magnetization will remain only if the domain to which the 
polarized droplet belongs is not eliminated. Then the re- 
manent magnetization will decay by the further domain 
growth as, 



SM n {t) ~ 5M n (t+) 2P t 



for 



'p(O) 



t > tf, 



T(L n /L )' 



rC 



- 1 



(81) 



Here C(L(t)/L) is related to the probability P s (L(t)/L) 
that a given site belongs to the same domain at the two 
different epochs L n and L(t)(> L n ) as given in (|73|). This 
is the aging part of the response which violates the FDT 
@ strongly. 

Combining the above results we can now evaluate the 
ZFC and TRM susceptibilities defined in (0) and (§) by 
summing over the responses at different epochs given by 
(|79|) and d8l]). Suppose that the observation time t be- 
longs to the n-th epoch and the waiting time f w belongs 
to the m-th epoch such that m < n. Then the ZFC 
susceptibility is obtained as, 



Xzfc(Mw) 



/ dt'R(t,t') = -J2SM k (t) 

k=m 



C [^](82) 



~ X (L(r),L(t)) 

„ f L(t) dL p{0)m 2 
+C i L(tw) L T(L/L y 

The first term in the last equation is the response of 
droplets (|79| ) during the last epoch where the observa- 
tion is being done and second term is the aging part ( |8l| ) 
due to the remanence of the response made at previous 
epochs. Note that the expression is valid also for the 
quasi-equilibrium and crossover regimes because there 
L(t) ~ L(t w ) holds and the second term is absent so 
that the expression becomes the same as (|54|). 

In the special case of t w = 0, the above result becomes 



XzFc(i,0) 



x(L(t),L(t)) 

" L{t) dL p{0)m 2 



^(Lit) 



rC 



XD 



L T(L/L Q ) e V L 
p(0)m 2 



where 



with 



T(L(t)/L y 



A{1) - C nst 



Cnst = c / dyy 
h 



-l-e+x 



> 0. 



(83) 



(84) 



(85) 



To obtain the last equation, we used the dynamical 
susceptibility \~d defined in mu) and the scaling form 
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C{x) ~ x~~ given in ( |75| ) and assumed that i is large 
enough so that L(t)/L -C 1 holds. Note that c nst repre- 
sents the strength of the aging part of the response whose 
contribution to xzFc(t,0) scales as c ns tm 2 /(L(t)/L ) e ■ 



Now comparing the result (|83|) with ( |10| ) which reads as 
Xfc = lim Xzfc(*,0) 

t — >00 

we find the dynamical susceptibility xd defined in j60| ) 
is nothing but the desired FC susceptibility xfc- Thus 
using the inequality (pi) we obtain, 



Xfc = Xd > Xea 



(86) 



which is nothing but the anticipated inequality (fl6|). 

Lastly let us evaluate the TRM susceptibility in the 
aging regime. Suppose the waiting time t w belongs to the 
m-th epoch and observation is done well after the waiting 
time such that the domain growth proceeds appreciably 
L{t) > L{tw). Then summing over the aging part of the 
response (J8l|) we obtain, 



XTRM^ = T + t w ,t w ) 



dt'R(t,t') ~^5M k {t) 



fe=i 



[ L ^dL ~p{Q)m\~(L{n 



L ?(L/L ) 

m 2 ( Lit) 

; T(L(t w )/i )« \L^j 
for Lit) ~ L(t) > Lit w ). 



^nst " 



L 

— A 



(87) 



Here c ns t is defined above in ( |S5| ) which represents the 
strength of the integral of the aging part of the response. 
To obtain the last equation, we used the scaling form 
C(x) ~ x~ x given in ([75]) and assumed that i w is large 
enough so that Lo/L(t w ) -C 1 holds. The functional forai 
agrees with what was anticipated by Fisher and HusetLa. 

Thanks to the sum rule (0), the relaxation of TRM 
susceptibility (||) and ZFC susceptibility (0) can be ob- 
tained from each other using xzfc(^ 0) obtained as (83). 
Here one must pay attention to the fact that xzFc(i,0) 
contains both the response due to droplets (which satisfy 
the FDT) and aging part (which violates the FDT) as 
can be seen in (B3). 



In the previous subsections IV E 2 and IV E 3, we ob 



taincd scaling properties of the ZFC susceptibility in the 
quasi-equilibrium/crossover regime. The ZFC suscepti- 
bility in these regimes contain only response of droplets 
which satisfy the FDT. However one can check that the 
TRM susceptibility in these regimes, which can be readily 
obtained via the sum rule, becomes a mixture of response 
due to droplets and aging part. Thus the ZFC suscepti- 
bility in the quasi-equilibrium/crossover regimes is better 
suited to examine responses of droplets than the TRM 
susceptibility. 

Conversely the TRM susceptibility in the aging regime 
is better suited to examine aging part of the response 
than the ZFC susceptibility. As can be seen in (|87j), the 



TRM susceptibility in the aging regime contains only the 
aging part of the response. On the other hand one can 
check that the ZFC susceptibility which can be readily 
obtained by the sum rule becomes a mixture of response 
due to droplets and aging part. 



6. Summary 

Let us consider large time limit such that L(t w ) 
is taken with the ratio 

x m 

Lit w ) 

being fixed to certain values. Here it is useful to consider 
the large time limit of the sum rule (|^). For any r that 
may be allowed to grow with t w we find, 

lim [xtrm(t + t w , t w ) + xzfcx(t + t w , tw)} = Xd 

(89) 

In the last equation we used the definition of the field 
cooled (FC) susceptibility defined in ( |l0| ) and our result 
(p6|). The asymptotic behaviors discussed below are dis- 
played in Fig. H and Fig. ||. 

First x < 1 corresponds to the quasi- equilibrium regime 
where the spin autocorrelation function slowly decays 
from 1 to static order parameter qfa accompanying some 
weak waiting time dependence or weak violation of time 
translational invariance (TTI). Here if one fix L(t) and 
let i(i w ) - * oo one obtains the ideal equilibrium limit 
behavior where the weak violation of TTI is removed. 
On the other hand, the FDT (|l|) 1 - C(r + t w ,t w ) = 
k^TxzFci T + ^wj^w) is satisfied even in the presence 
of the weak violation of TTI. In the large time limit 
Lit w ) — > oo with fixed x < 1, the time dependences (in- 
cluding the weak violation of TTI) disappears such that 
the spin autocorrelation converges to the static EA order 
parameter q^A and the ZFC susceptibility converges to 
the equilibrium susceptibility k^TxEA- Correspondingly 
the TRM susceptibility converges to k^Txu — ^b^Xea = 
<7d — 5ea because of the sum rule ( |S9| ) . In the last equation 
we used k^Txr> = 1 — Qx> given in (60). To summarize, 
the spin autocorrelation function and the susceptibilities 
asymptotically become flat lines in the quasi-equilibrium 
regime (x < 1) as displayed in Fig. ^| In the paramet- 
ric plot of Fig. [|, the whole quasi-equilibrium regime 
converges to a single point (qfaj ^b^xea)- Importantly, 
the scaling theory has provided not only such asymptotic 
limits but details of finite-time (r, t w ) correction terms by 
which t he asy mptotic limit is approached as discussed in 
section IV E 2 . The latter are well amenable to be exam- 
ined seriously by experiments and simulations. In section 
VII we present such a detailed analysis for the 4D EA 
model based on MC simulations. 

Second x ~ 1 corresponds to the crossover regime. 
From (|69[), we expect the spin autocorrelation function 
drops vertically (against L(r)/L(t w )) from the EA order 
parameter ^ea down to the dynamical order parameter 
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5d which is smaller than <7ea due to the anomalously soft 
droplets. In this regime we expect the FDT ( |l8| ) is still 
satisfied in spite of the strong non-stationary character 
and we expect the ZFC susceptibility increases vertically 
from the static susceptibility k^TxEA to a larger value 
kBTxB = 1 — (ZD- Correspondingly the TRM susceptibil- 
ity drops off vertically from k^TxD — k B TxEA = Qea — <Zd 
to zero because of the sum rule (^9|). In the parametric 
plot of Fig. ||, the crossover regime converges to the 
line of points in the section of the FDT line between 
(<?ea,&bTxea) and (? d ,A;bTx d ). 

The abruptness of the changes of the two time quan- 
tities at x ~ 1 is very surprising. Here let, us recall the 
well known experimental^ and numericalE-3 observations 
that the so called relaxation rate 



S(t, t w ) 



_ cIxzfc(t + t w ,t v 
dln(r) 



(90) 



has a broad peak centered at around t ~ i w . Our sce- 
nario naturally suggests a modified relaxation rate func- 
tion, 

S mo d{x,t V! )= . (91) 

X x=L{T)/L(t„) 

which should have a sharper peak at x ~ 1 with increas- 
ing L(t w ). This modified relaxation rate will be useful 
for further numerical simulations and experiments. 

In order to describe the interior of the crossover regime 
more closely, different scaling variables other than L(t) 
and L(t w ) are certainly needed. Unfortunately, the 
droplet theory which is based on the dynamical length 
scales cannot provide information for proper scaling vari- 
able to describe the interior of the crossover regime. A 
possible scaling variable would be a = r/t w as proposed 
by Fisher and HuseO-3 with which the abruptness will be 
absent. Note that possible large time limits t w — > oo 
classified by different values of a are smashed to a point 
x = L(t) I L(t w ) = 1 in the x-axis for all a. 

Finally x > 1 corresponds to the aging regime where 
the spin autocorrelation function takes a continuous 
value between go and as deceasing function of x as 
given in @. Here the result @ implies the TRM 
susceptibility is asymptotically zero. The latter means 
the ZFC susceptibility converges to the FC susceptibility 
k B TxFC because of the sum rule @. The FDT @ is 
thus strongly violated. In the parametric plot of Fig. [|, 
asymptotic limit for each x converges to a point on the 
flat horizontal line connecting (0, Ixd) and (qu, ^b^Xd)- 
In conventional u ndersta nding)!! which includes usual 
3 ffl and 



ing systems!: 



mean-field spin-glass 



coarser. 

modelso, both the violation of TTI and FDT happens 
asymptotically at the same static oder parameter ^ea- 
A remarkable feature of our present scenario is that the 
break points of TTI and FDT are separated: the break 
point of TTI is located at q^A while that of FDT is lo- 
cated at the dynamic oder parameter . 



So far it is implicitly assumed that FDT (18) is valid 
in the quasi-equilibrium regime and also in the crossover 
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FIG. 2: Different asymptotic regimes of the two-time quanti- 
ties. See text for details. 
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FIG. 3: Different asymptotic regimes in the C — T\zfc plane. 
The dotted line labeled k B Tx = 1 - C represents the FDT 
line. See text for details. 



regime in spite of the fact that there are weak and strong 
waiting time effect. A supplementary argument for the 
validity of FDT can be made by considering the bound 
on the possible violation of FDT found inE3. The rigor- 
ous boundE3 on the integral violation I(t,t w ) defined in 
( p3| ) is put in terms of the entropy production rate which 
presumably has the same scaling form as the energy re- 
laxation rate. Then using the scaling form of the energy 
relaxation (B7f) one finds, 



\I(t,U)\ < K J d S] /^ 

:> ((t/i ) 1/2+Q 



< K' 



(W*o) 



l/2+c 



(92) 



with K and K' being certain finite constants and < 
a <C 1 being an arbitrary small positive non-zero num- 
ber. Now let us consider the bound in the large time 
limit t w — > oo in the quasi-equilibrium regime x = 
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L(t) I L(t w ) < 1. For convenience let us introduce y such 
that, 



t — t w — to(t w /to) V 



(93) 



Then we find 



lim I(t,t w ) < (1/2 + a)K' lim (t w /t Q )- 1/2+v+a 



= 



(94) 



The last equations holds for < y < 1/2 since a is an ar- 
bitrary small positive number. This observation implies 
that FDT is satisfied not only in the equilibrium limit 
but also in the quasi-equilibrium regime. In order to ver- 
ify the validity of FDT in the crossover regime y ~ 1, 
apparently improved bounds are needed. 



7. Comparison with Conventional Pictures 

Finally, let us compare our scenario with the conven- 
tional picture for isothermal agingQ which applies for the 
dynamical mean-field theoigi-|tMFT) of spin-glassesEl and 
usual coarsening systemsOElMiil. Here let us consider the 
asymptotic limit t w — > oo of the two time quantities by 
fixing the value of the autocorrelation function to a cer- 
tain value C in the range < C < 1. In that limit, the 
ZFC susceptibility becomes a function of C, \ = x(C)- 
This definition of asymptotic limit can be considered in 
general and our scenario implies the simple structure of 
x(C) shown in Fig. |^. In the usual coarsening systems 
and in the dynamical MFT, the FDT line terminates 
at (<7ea,Txea) wn ^ e ^ extends up to (<Zd,Txd) in our 
scenario. In the usual coarsening systems, a horizontal 
line connects (0, /cb^Xea) and (<?ea, ^b^Xea) where the 
FDT is violated strongly. On the other hand, the dy- 
namical MFT predicts a curved line between (0, xd) and 
(<7ea> Xea) ■ Our picture is different from both of them. 

The difference between usual coarsening systems and 
our scenario is the presence of the anomalously soft 
droplets. Droplet excitations .exist also in simple coarsen- 
ing systems like ferromagnetsE-l but with extremely small 
probability. In our scenario, the anomalously soft droplet 
can exist with probability of order O(l) in a given domain 
irrespective of the size of the domain so that the ther- 
mal fluctuations in each domain are anomalously large 
compared with equilibrium where there are no extended 
defects. _ 

For clarity, let us notecj that there are also some 
excessive response in usual coarsening systems due to 
thermalized domain wallscil at wavenumbers k such that 
kL(t) ^> 1. They are similar to the anomalously soft 
droplets in our scenario in the sense that 1) they satisfy 
FDT but 2) disappear in the ideal equilibrium. However 
their integral contribution to the response decreases with 
the growth of the domain L(t) so that their contribution 
vanishes asymptotically. On the other hand the excessive 
response of the anomalously soft droplets in our scenario 



give O(l) contribution irrespective of the size of the do- 
main so that their contribution to the response do not 
disappear as far as t w — > oo limit (ideal equilrbium) is 
not took first. 

Although both the dynamical MFT and our scenario 
concludes the inequality ( |l6|) xfc > Xea, the origins 
are different. The difference between the two, called 
anomaly, is attributed to contributions of aging part of 
the response function which violates the FDT in the case 
of the dynamical MFT while it is rather attributed to 
the responses due to the anomalously soft droplets which 
still satisfy the FDT in our scenario. For example, one 
can see in ( |83] ) that xfc(= Xd) is generated not from the 
aging part but from the part due to the anomalously soft 
droplets which satisfies ihe FDT. 

There is a conjectureE^I that there is a connection be- 
tween the statics and dynamics such that the function 
X = x(C) i n the dynamics should be related to the over- 
lap distribution function P(q) in equilibrium, as 



fcsTx(C) 



dC I d,C"P(C"). 



(95) 



This is -known to hold exactly in some (not -all) mean-field 
modelsa and usual coarsening systemsESEJ. Moreover, 
an interesting conjecturalj was proposed recently that at 
finite time (length) scales, the above formula holds for 
the two time quantities at a finite waiting time t w and 
a P(q) measured in equilibrium of a finite system of size 
L(t w ). We agree with this proposal partially but not 
completely as we discuss in the following. 

In order to fully reproduce the parametric plot Fig. 
|3|, we would need P(q) with a delta peak at q = q^ 
rather than at q = qEA- On the other hand, we expect 
a P(q) with delta peak at q — <?ea in the absence of 
any extended defects as the original droplet theory has 
predicted!-! Thus the correspondence between the dy- 
namics and statics does not hold at all in this sense in 
our scenario. 

However, one can explicitly consider a rather special 
static situation in the presence of extended defects of 
size, say R, which is actually what we considered in sec- 
tion IV D. In any finite size systems, it is likely that 
the existence of boundaries (periodic/free e.t.c.) will in- 
trinsically ind uce . pertain defects as compared with in- 
finite systemsE3c3'c3. Thus we expect actually the cir- 
cumstance we are considering is relevant in practice. Our 
scenario implies the average overlap q = J dqqP(q) mea- 
sured in such an equilibrium is equivalent to q(R, R) = 
q D + p(0)m 2 k B T/T(R/L ) e A{l) of (|f). Thus our sce- 
nario suggests xfc = lim T _,aoXzFc(T + i w >*w) = l-q = 
l — f dqqP(q) for any t w while linv^oo C(r+t w , t w ) — > 
for any t w . Thus our scenario agrees with the conjecture 
of Ref. [l5] at the special point C — plus the usual FDT 
part between C = <7ea and C = 1 but not in the sec- 
tion < C < qEA- Here P(q) should be understood as 
measured in equilibrium but with the extended defects. 
It will be natural to expect that P(q) in such a situ- 
ation becomes a non-trivial, non-self-averaging function 
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of q with finite amplitude at q — as obsepjed. in many 
numerical simulations of finite size systemsHH 

In the above arguments, we considered asymptotic 
limits with fixed C while it is more natural to con- 
sider asymptotic limits with fixed ratio of the two length 
x = L(T)/L(t w ) in the droplet theory. In the dy- 
namical MFT, one needs infinitely many kinds of time 
reparametrization functions h(t) to span the whole cor- 
relation range < C < 1 but the droplet theory has 
only one natural variable L(t). However it should be re- 
marked that the scaling variable to describe the interior 
of the crossover regime qu < C < q~EA is not known at 
present within the droplet theory. 



GROWTH LAW OF THE CORRELATION 
LENGTH IN OFF-EQUILIBRIUM 



As discussed in Sec. IV, the dynamical length scale Lit) 
and its growth law with time are important to understand 
nature of aging in SG systems from the view point of 
the scaling theory. In the present section, we discuss 
our numerical results of the length scale and its growth 
law during isothermal aging. A plausible definition of 
the length scale L(t) is given by a decay constant of a 
correlation function G{r,t). We measure the equal-time 
replica correlation function in off-equilibrium under zero 
magnetic field, defined by 

G{r,t) = ^(Sl a \t)sW(t)S$(t)sW(t)), (96) 

i 

where a and f3 denote the replica indices which are up- 
dated independently from different initial random spin 
configurations. In our simulation, only one MC se- 
quence is performed for each random bond. Typical 
number of samples averaged over bojad-jrealizations is 
about 128. As reported previouslyErE3E3, the correla- 
tion function G(r, t) exhibits a complicated functional 
form which is not a simple exponential and depends on 
time t. This may be because a characteristic distance 
above which G(r, t) approximately follows an exponen- 
tial form depends on time considerably. In order to avoid 
the artificial effect, we do not employ the so-called ratio 
methocOO where the full data of G(r, t) are used inde- 
pendent oft. Instead, we estimate L(t) by fitting directly 
the tail part of G(r, t) to an exponential formula for each 
time t. In the fitting procedure, we focus our attention 
only on the large distance tail of G(r, t) and carefully 
choose the range depending on time t. 

In Fig. [| we show our results of the length scale L(t) 
at the critical temperature T c = 2.0J. The data of the 
length scale L(t) follows expected critical power law ( |30| ) 
except for the data with L = 8, which is due to finite-size 
effect. The estimate of the exponent z, 4.28(5), is roughly 
consistent with that of the previous workE3, which quoted 
z = 4.45(1). The data, including the size L = 8 deviated 
from the power law, scale well using a standard finite-size 
scaling form, shown in the inset of Fig. [l| It is seen that 
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FIG. 4: R(t) of the Ad Ising SG model at T/J = 2.0 with 
different system sizes. The dashed line represents a power 
law with exponent 1/z. In the inset, an expected finite-size 
scaling plot is shown. 
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FIG. 5: Time evolution of R(t) of the Ad Ising SG model at 
T c and below T c with L = 24. 



the length scale Lit) manages to reach at most few lattice 
spacing even at about 10 5 MCS. Nonetheless it should be 
noted that L(i) already captures macroscopic behavior 
in the sense that the critical exponent z is successfully 
estimated from Lit) and that a strong finite-size effect is 
observed already in L = 8 data. 

Let us now examine the crossover from critical to ac- 
tivated dynamics discussed in section IV Bj . In Fig. ||, we 
display the raw data of the length scale Lit) below T c and 
at T c . The critical dynamics is expected to be dominant 
near T c . In fact, the length scale Lit) at T/J = 1.8 is 
not distinguishable from that at T = T c within our time 
regime. Let us examine the crossover scaling defined in 
( p4"| ) and (|35"|). We present the scaling plot in Fig. ^by us- 
ing the same data as shown in Fig. || In this plot, for two 
scaling parameters,!-^ and v, we use the known values 
obtained previouslyc^l, while only one remaining param- 
eter tjj is appropriately chosen for the data with different 
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temperatures to merge into a universal curve. The best 
scaling plot is obtained by tp ~ 2.5 — 3.0. The proposed 
scaling pretty well works in the observed time regime. 
It is clearly found that the scaling function exhibits a 
crossover from the critical power law to slower growth 
law associated with low temperature dynamics which is 
compatible with the logarithmic growth law predicted by 
the droplet theory. 
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FIG. 6: Scaling plot of R(t) of the 4d Ising SG model, where 
the SG transition temperature T c = 2.0 and the critical ex- 
ponent v — 0.93 are fixed, but the dynamical exponent z is 
obtained to be 4.98(5) from the best scaling. 



VI. RELAXATION OF ENERGY AND DENSITY 
OF DOMAIN WALL 

In the present section, we examine scaling properties of 
one-time quantities under isothermal aging after quench 
based on the view point presented in section IV C . In the 
following sections, we perform simulations mainly at two 
low temperatures T/J = 1.2 and 0.8 which amount to 
0.6T C and 0.4T C , respectively. It is found that the effects 
of critical fluctuations do not dominate our time window 
at these temperatures but can be taken into account in 
a renormalized way. 



A. Energy 

In Fig. 0, we present the data of the energy per spin 
e(t) defined in (||) at T/J = 0.8 and 1.2. As seen in 
the upper figure, the energy function relaxes with time 
to an equilibrium value at each temperature. However, 
it is rather hard to extract the equilibrium value from 
the figure because of the extremely slow dynamics. On 
the other hand, the scaling formula (|37]) says that the en- 
ergy approaches its equilibrium value linearly as a func- 
tion of L (t). We thus plot the same data set against 
L e ~ d (t) in the lower one by using the length scale L(t) 
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FIG. 7: Energy relaxation at T / J = 0.8 and 1.2 as a function 
of elapse time (upper figure) and the length scale L(t) (lower 
figure). The straight lines in the lower figure represent fitting 
result to a linear function function of (L(t)/lo) 6 ~ d . 



estimated independently in the previous section and the 
known value of the stiffness exponent 6 = 0.82E3. We 
see the linear behavior as a function of L e ~ d (t) for large 
L(t) which supports the validity of the scaling formula 
(p7|). One can take the long time limit of the energy 
relaxation using the scaling formula. A similar scaling 
analysis of the energy has been confirmed ip_the 3d Ising 
EA SG model including finite size effectstj. This two 
stroke strategy has an advantage over the direct analy- 
sis with time and gives us a more powerful tool in the 
analysis of two-time quantities discussed in the following 
sections. 



B. Density of domain wall 

In FigJS, we show the domain- wall density p s (t) de- 
fined in (B9). The average over bond realizations is taken 
over 256 samples with L = 24. It is found that the den- 
sity monotonically decays and there is no tendency of 
saturation, which, is similar behavior observed in the 3d 
Ising EA modelEil. This is clearly seen in the lower fig- 
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FIG. 8: Density of the domain wall at T/J — 0.8 and 1.2 with 
L — 24 as a function of time f(upper figure) and L(t) (lower 
figure). The straight line in the lower figure representp-fitting 
result to a linear function of L d f~ d (t) with d f = 3.75E1 



ure, where following ( p8[ ) we plot p s {t) as a function of 
L d f~ d (t) estimated in the previous section. Here we use 
the value of the fractal dimension df = 3.75 recently eval- 
uated in Ref. As shown in the lower figure, p s (t) is 
well fitted by a linear function of L d f~ d (t) and the fitting 
function is down to zero in the large time/length limit. 
We stress again the validity of the scaling formula with 
the length scale L(t). 



VII. RELAXATION OF CORRELATION 
FUNCTION AND LINEAR SUSCEPTIBILITIES 

We now turn to the two time q uantiti es. First we dis- 
play the data in the next section VII A and d iscuss the 
overall features qualitatively in section VII B from the 



point of view of the scaling theory explained in section 
EVE. Then we go to more detailed examinations of the 
scaling properties in the subsequent sections. In order 



to test the scaling ansatz explained in section [V E which 
are expressed in terms of the dynamical length scale L(t), 
we use the data of L(t) discussed in the previous section 
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FIG. 9: Spin autocorrelation function of the 4D Ising SG 
model at T/J = 1.2 (upper figure) and 0.8 (lower figure) with 
different waiting time. The top data curves are the equilib- 



rium curves obtained in section VII C 



A. Measurements of Spin Autocorrelation 
Function and Linear Susceptibilities 

In Fig. ||, we present the data of the spin autocorre- 
lation function C(r + t w ,t w ) measured up to 10 5 MCS 
for various waiting times t w = 10, 10 2 , 10 3 , 10 4 , 10 5 at 
T/J = 1.2 and 0.8 . The system size is L = 24. The data 
arc obtained by performing MC simulations starting from 
random initial conditions. The average over realizations 
of randomness is taken over 32 samples. The data show 
clear waiting time dependences. We also show the curve 
of equilibrium relaxation C oq (r) extrac ted in the analysis 
which will be explained later in section VII C . The latter 
yields the value of static EA order parameter (?ea whose 
value is also indicated in the figure. In addition we also 
indicate in the figures the values of the dynamical or- 
der parameter go which will be obtained in the analysis 
explained later in section VII E. As discussed in section 



IV E, we expect that the spin autocorrelation function de- 



velops a plateau at gEA in the quasi-equilibrium regime 
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and decays down to zero in the aging regime. The decay 
is expected to be most steep at around the dynamical or- 
der parameter <7d • The data indeed appears qualitatively 
compatible with these expectations. 

In Figs.|H^ and O, we display the data of the ZFC 
susceptibility (T/J]xzfc( t + t w ,t w ) and the TRM sus- 
ceptibility (T I J)xtrm(t + tw, tw) measured up to 10 5 
MCS for various waiting times i w = 10, 10 2 , 10 3 , 10 4 , 10 5 
at T/J = 1.2 and T/J = 0.8 using field of strength 
h/J = 0.1. Again the system size is L = 24. The average 
over realizations of randomness is taken over 32 samples. 
Here the susceptibilities are defined by dividing the mea- 
sured magnetization per spin m(t) = (1/N) 5^2=1 &i (*) 
at time t by the strength of the external magnetic field 
h/J. For the measurement of the TRM susceptibility 
Xtrm (t = t + t w ,t w ), some waiting time t w is elapsed 
under the field and then the field is switched off at t = t w . 
Conversely, the ZFC susceptibility xzFc{t = t + t w ,t w ) 
is measured by first elapsing the waiting time t w with no 
applied field and then field is switched on at t = t w . 

To check the linearity of the response, we examined the 
sum rule (^|) : the sum of the two susceptibilities become 
(T/ J)xzFc(t, 0). As shown in FigjlJ, the sum rule is 
satisfied over our time window within the statistical ac- 
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FIG. 11: Linear susceptibilities at T/J — 0.8. Here 
(T/J) X ea = 0.2 and (T/J) X o = 0.33. 
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FIG. 12: Check of the sum rule (|) at T/J = 1.2 (upper fig- 
ure) and 0.8 (lower figure). The data of the ZFC susceptibility 
{T/J)xzFc(t, 0) is also included. 



curacy. 

In Figsjl?] and O, we can see clear waiting time de- 
pendences of the susceptibilities. As discussed in section 
IV E, we expect that the ZFC susceptibility first devel- 



ops a plateau at the static susceptibility xea within the 
quasi-equilibrium regime and then grows further up to 



19 



the dynamical susceptibility xd later in the aging regime. 
Correspondingly, we expect that the TRM susceptibil- 
ity develops a plateau at xd ~ Xea within the quasi- 
equilibrium regime and decays down to zero in the aging 
regime. For the references, we indicated the values of the 
static susceptibility (T/J)xea , the dynamical suscepti- 
bility (T/J)xd and their difference in the fig ures u sing 
the va lues which will be obtained in section VII C and 
VII E] , The data indeed appears qualitatively compatible 
with the expected behavior. 



B. Overall Features 

In Fig. [D|, we show the spin autocorrelation function 
C(t + f w , i w ) and the ZFC susceptibility 1 - (T/ J)%(t + 
t w ,t w ) plotted against x = L(t) / L(t w ) . 

Here we parametrized the times r and t w using the 
time dependent length L(t) obtained in section [v|. More 
precisely, we first fitted the data of the domain size L(t) 
as L(t)/L = a Q + ax log(i) + a 2 log 2 (t) + a 3 log 2 (i) with 
t measured in the unit of MCS and Lq = 1 (lattice dis- 
tance). This fitting was enough to model the data of L(t) 
within our time window. Then we used the resultant fit- 
ting functions to do the parameterization here and in the 
following analysis. When we do the parameterization, 
we discard short time data t < 10 because L(t) is not 
available there. 

The figure should be compared with Fig. || which ex- 
plains the expected asymptotic behavior of the two time 
quantities in the large time limit L(t w ) — > oo with the 
ratio x = L(r)/£(i w ) fixed to certain values. We ex- 
pect three distinct regimes depending on the value of 
the fixed ratio i) quasi-equilibrium regime x < 1 where 
the autocorrelation functions spans values in the range 
<7ea < C < 1 ii) crossover regime x ~ 1 for qn < C < <7ea 
and iii) aging regime x > 1 for < C < qr>- Here q&A 
and gD are convergence points of the break points of TTI 
and FDT respectively in the limit L(t w ) — * oo. 

In the quasi-equilibrium regime x < 1, we expect 
that the spin autocorrelation function convergence to 
a plateau <7ea in the large time limit L(t w ) — > oo. 
Correspondingly the ZFC susceptibility (T/J)xzfc( t + 
t w ,t w ) is expected to converge to the static susceptibility 
(T/J)xea = 1 — Qea- In the figure, we indicat ed the value 
of (?ea which will be obtained later in section VII C. The 



data indeed appears descending down toward the plateau 
from above with increasing £ w (at x < 1) but still far from 
it. Fortunately, the scaling theory provides prediction on 
the correc tion ter ms to the asym ptotic limit as explained 
in section [VE2. In the section VII C, we will examine 



the correction terms in detail which actually yields the 
value of <7ea- Second important observation is that the 
FDT 1 - C(Mw) = (T/J)xzfc(Mw) is well satisfied in 
the quasi-equilibrium regime x < 1 as expected. 

In the crossover regime x *~ 1, we expect a vertical 
drop from the plateau at the static EA order parameter 
<7ea down to the dynamical oder parameter for the 
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FIG. 13: Plot of spin autocorrelation function C(r + i w ,t w ) 
and 1 — {T/J)xzfc(t + iw,iw) against x = L(r)/L(i w ) at 
T/J — 1.2 (upper figure) and 0.8 (lower figure). For each 
waiting time t w , the lower curve is the spin autocorrelation 
function and the upper curve is the susceptibility. 



spin autocorrelation function. Correspondingly the ZFC 
susceptibility [T / J)xzfc( t + iw,^w) is expected to jump 
from the static susceptibility (T/J)xea = 1 — <?ea up 
to the dynamical susceptibility (T/J)xu = 1 — qu still 
keeping the FDT 1 - C( t, t w ) = (T/J)xzfc(*,*w) as ex- 
plained in section IV E 2} In the figure, we indicated the 
value of which will be obtained later in section VII E . 
We will exam ine the slope of the curves around x ~ 1 in 
section VII E and find indeed that the suitably refined re- 
'mod(x,t w ) defined in ( |9l| ) have a 
1. Furthermore, 



laxation rate function S, 
sharply pronounced peak at around x 



we will examine the violation of FDT I(t + i w ,£w) 



(T/J)xzFc(r + i w ,tw) defined in (E 



l-C(r + t 

in section VII G and find indeed that it is decreasing with 
increasing f w in the crossover regime x ~ 1. The result is 
compatible with our expectation that the FDT is asymp- 
totically valid even in the crossover regime. 

In the aging regime x > 1, we expect the ZFC sus- 
ceptibility 1 — {T/J)xzfc(t + t w ,t w ) converges to the 
dynamical susceptibility (T / J)xd = 1 — <tD while spin 
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FIG. 14: Parametric plot (T/J)xzFc(t,*w) vs C(t,t v ) at 
T/J = 1.2 (upper figure) and 0.8 (lower figure). The straight 
tangent lines represents the FDT ([l8]). The convergence 
points of the break points of TTI (qea, (T / J)xea) and the 
break points o f FDT f an, (T/J ) Yn) which will be obtained 
later in section VII C and VII E respectively are indicated. 



autocorrelation function C(t + t w , t w ) becomes a scaling 
function of L(r)/L(t w ). Indeed the data appears slow 
converges to such limits. For the susceptibility, the scal- 
ing theory provides prediction on the scaling forms of the 
correct ion term s to the asymptotic limit as explained in 
section ]FVE5| . We will examine and confirm them in sec- 
tion VII E . The value of qu indicated in Fig. [l^ will be 
actually obtained as the result of such an analysis. 

In Fig. the susceptibility (T/J)xzfc(t + t w ,t w ) 
and the spin autocorrelation function C(r + t w ,t w ) are 
plotted in a parametric way. It should be compared with 
Fig. which explains the expected asymptotic regimes 
in the large time limit L(t w ) — > oo with the fixed ra- 
tio x = L(r)/L(t w ). The curves apparently continue 
to move upwards with increasing waiting time t w and 
the break point of FDT does not appear to conv erge to 
(c7ea, (T/ J)xea) which will be obtained in section VII C . 
This observation is consistent with our picture presented 
in Fig. It is radically different from conventional 
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FIG. 15: Spin autocorrelation function at T/J — 1.2. The 
top data points are the equilibrium limit obtained by the ex- 
trapolation method explained in the text. 
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FIG. 16: Typical examples of the extrapolation to the long 
t w limit with fixed r. Two arrows represent the extrapolated 
values of C e ,(r) of r = 100 and 1000 at T/J = 1.2. The 
horizontal error bars are those of L(t) presented in Fig. H. 



understandingu that assumes break points of TTI and 
FDT are the same. As we discuss in section VII F, the 



separation of the break points of TTI and FDT will be- 
come apparent when time (length) dependences are ex- 
plicitly considered. 



C. Quasi-Equilibrium Regime 

Let us now examine the scaling properties in the quasi- 
equi librium regime x = L(t) / L(t w ) < 1 discussed in see- 
Since we have seen the FDT is well satisfied 



IVE2 



tion 

in the quasi-equilibrium regime x < 1 (see Fig. |l^), we 
will only use the data of the spin autocorrelation func- 
tion. To focus on the quasi-equilibrium regime, we took 
another dense data set of the autocorrelation function 
C(r + t w ,t w ) (shown in Fig. 15) of many waiting times 
< w = 1, . . . , 30000(MCS) and relatively short time sepa- 
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FIG. 17: Equilibrium spin autocorrelation function C e q(r) at 
T/J = 1.2 (upper data)and 0.8 (lower data). The horizontal 
error bars are those of L(t) presented in Fig. ^. The lines 
represent fittings according to the formula (pq). 
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FIG. 18: Scaling plot of the non-equilibrium correction terms 
of C(t + t v ,t m ) in the quasi-equilibrium regime at T/J — 1.2 
(upper figure) and 0.8 (lower figure). The horizontal error 
bars are those of L(t) presented in Fig. 0. The line has the 
expected slope d — 9 for small L(t) / L(t lL !) limit (see (^)). 



rations r < 1000 (MCS). Here we don't use the fitting of 
the growth law of L(t) to parameterize the time but use 
directly the data of L(t) obtained in section [v|. 

The formula ( |67| ) combined with (|6^) gives a protocol 
to extract the static EA order parameter as the following. 
This ..analysis is already reported partly in our previous 
workED. 

At first we take the equilibrium limit t w — > oo of C(r + 
t w ,t w ) for each time separation r. From (|67|), one finds 
that the autocorrelation function becomes only a linear 
function of 1 / L d ~ e (t w ) for a fixed r. Thus we plotted the 
data points of a given r of various t w against 1/L d ~ e (t w ) 
and fitted to a linear function in the large L(t w ) regime. 
The equilibrium limit C cq (r) is read off directly ftop the 
fit. Here we use the stiffness exponent 8 = 0.82c3. For 
L(t w ) we used the data obtained in the previous section 
Typical fitting results are shown in Fig. |l^. The 
linearity as a function of 1/L d ~ e (t w ) supports the validity 
of the formula (|67|). We repeated the same procedure 
for each r. The obtained equilibrium curve C eq (r) was 
displayed in Fig. |[ 

Next we take the large r limit of the extracted C eq ("r) 
using (poj). As shown in Fig. |l7], C cq (t) appears as a 
linear function of 1 / L e (r) which support (|65|) . A linear fit 
gives the EA order parameter <7ea as the limiting value. 
We obtained the EA order parameter ^ea = 0.58(2) at 
T/ J = 1.2 and 0.82(1) at TUJ = 0.8. To our knowledge 
this and our previous workE-3 is the fist which confirmed 
this fundamental scaling law (|65|). 

Finally, we discuss the weak non-equilibrium correc- 
tion term to the equilibrium limit which is responsible 
for the weak violation of TTI in the quasi-equilibrium 
regime. The expression (^) suggests that the weak non- 
equilibrium correction term AC(r + t w ,t w ) — C(t + 
tw, tw) — Ceq(T) multiplied by L (t) becomes only a func- 
tion of L{r)/L(t w )(< 1). As shown in Fig. [T^, we confirm 
this scaling form. For the limit of L{t) / L(t w ) <C 1, the 
scaling function shows the expected power law behavior 
(L(t) I 'L(t w )) d ~ e being consistent with (|67 



D. Crossover Regime 

In the crossover regime x = L(t) / L(t w ) ~ 1, we 
expect a vertical jump of the two time quantities at 
asymptotic limit L(t w ) — > oo. In Fig. [l9] we display 
the plot of the usual relaxation rate function S(r,t w ) 
of (p0[), which shows the well known peak structure ob- 
served in experiments^] and a MC simulatiora. This 
can be interpreted as the signal of the rapid changes in 
the crossover regime. However, our scenario naturally 
suggested a more appropriate, modified relaxation rate 
function S mo d(x,t w ) of (pij). In Fig. gO|we show the plot 
of the modified relaxation rate function. It clearly devel- 
ops a sharp peak at around x ~ 1 with increasing L(t w ) 
as expected. 
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E. Scaling of Susceptibilities 

For the growth of the ZFC susceptibility with zero 
waiting time we expect a scaling form (83) which reads 



Xzfc0,0) ~ xd 



c"'m 2 



T(L(t)/L y 



with d" being a numerical constant defined in (p4j). 

In Fig. |l| we show the sum (T / J)xzfc(t + t w , t w ) + 
(T/J)xtrm(t -Mm,tw) plotted against L(t = r + t w ) _e 
using = 0.820. We have checked that they satisfy 
the sum rule @ and agrees with (T/J)xzfc(A0) (see 
Fig.|l2]) . As can be seen in the figure, the data are in- 
deed consistent with the expected scaling form being lin- 
ear with 1/L(t) and pointing toward a constant in the 
limit L(t) — > oo. From the linear fit shown in the fig- 
ure we find the values of the dynamical susceptibility 
(T/J)xd as 0.54 at T/J = 1.2 and 0.33 at T/J = 0.8. 
The corresponding dynamical oder parameter can be 
determined via FDT {T/J) X d = 1 - <7d © as 0.46 at 
T/J = 1.2 and 0.67 at T/J = 0.8. 

In the analysis of section VII C , we have obtained 



from which we readily find the equilibrium susceptibility 
(T/J)xea = 1 — Iea- Interestingly enough we find the 
data of (T/J)xd(£, 0) shown in the figure clearly goes 
over the static susceptibility (T/J)\ea and the antici- 
pated inequality xfc = Xd > Xea holds (see ( |36"| ) ) . This 
is one of the main results of the present numerical simu- 
lation. 

The sum rule (|9|) requires us to examine only either the 
TRM or ZFC susceptibility plus the growth of the ZFC 
susceptibility with zero waiting time xzfc {t-i 0) which was 
obtained above. We have already analyzed the ZFC sus- 
ceptibility in the quasi-equilibrium regime and crossover 
regime in the preceding sections. In the following, we ex- 
amine the TRM susceptibility in the aging regime which 
will complete our scaling analysis of the linear suscepti- 
bilities. 

Let us examine the scaling ansatz (^) which reads as, 



XTRm(t + tw,t w ) ~ C nst 



L(t) 

T(L(t w )/L ) e \L(t w ) 



the value of the equilibrium EA order parameter ^ea 



for the TRM decay in the aging regime L(t) / L(t w ) > 1. 
In Fig. 22 , we show the scaling plot of the the decay of the 
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FIG. 21: Growth of the sum (T/J)xzfc(t + t w ,*w) + 
(T/ J)xtrm(t + t w , iw) with The data is the same as those 
used in Fig. [l2. The fitting lines are due to (|S3|). 



TRM susceptibility xtrm using the stiffness exponent 
9 = 0.82O. The results indeed agrees very welL_with 
the prediction ( |S7|) anticipated by Fisher and HuseEj and 
gives the non-equilibrium exponent A ~ 3.5. The latter 
satisfies the bound (|79) with d = 4. To our knowledge 
this is the fist time that this fundamental scaling law is 
confirmed. 



F. Separation of the Break points of TTI and FDT 

The results of the analysis presented so far supports 
well the existence of the anticipated asymptotic regimes 
displayed in Fig. || and Fig. || which predicts separation 
of the breaking of TTI and FDT. In order to look more 
directly on the separation, we present in Fig 
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_ the ZFC 

linear susceptibility (T / J)xzFc(t, t w ) and spin autocor- 
relation function 1 — C(i, i w ) of various waiting times t w 
plotted against 1/L(t) 9 . 

In the figure, we included the equilibrium ZFC linear 
susceptibility {T / J)xccSj) for comparison. The latter 
was obtained by the FDT (T/J)x eq (r) = 1 - C* cq (r) us- 
ing the equilibrium spin autocor relatio n function C eq (T) 
extracted previously in section VII C] . Since C cq (r) is 
a line ar function of 1/L(t) as we confirmed in section 
VII C (see Fig. [l7]), (T/J)x cq (r) becomes a straight line 
in the plot pointing toward the equilibrium susceptibility 
(T/J)xea- The top curve is (T/J)xzfc( t , 0) which is at 
the other extreme: zero waiting time. It becomes also 
a stra ight line in the plot as we already saw in section 
VII E] (See Fig. ^lj) pointing toward the dynamical sus- 
ceptibility (T/ J)xd which is significantly larger than the 
equilibrium susceptibility (T/J)xt>- 

Note that the slope of the susceptibility in the zero 
waiting time limit is smaller than that of the equilib- 
rium limit. This can be explained as the following. From 
(Bra) , the slope of the equilibrium limit is expected to be 
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FIG. 22: Scaling plot of (T/J)xtrm(t + i w ,*w) at T = 1.2 
(upper figure) and T = 0.8 (lower figure). The dotted line is 
a power law fit with exponent A = 3.5. 



proportional to c defined in ( p4[ ) while ( p3| ) suggests the 
slope of the zero waiting time limit is proportional to 
c'" = A(l) — c nst given in(p4]). Here c ns t is a positive 
constant and the value of A(l) can be evaluated by ( p?|) . 
The latter implies A(l) < c because the first term in ( |57| ) 
becomes identical to c defined in (^4j) while the second 
term is negative because T/T cf j[y] > 1. Thus we expect 
c > c'". 

Now a surprising observation is that the break points 
of FDT moves further away from the equilibrium curve 
by increasing t w . It appears very unlikely that the break 
points of the FDT converge to the equilibrium suscep- 
tibility (T/J)xea- This feature strongly suggests that 
the violation of TTI and FDT do not take place simul- 
taneously but separ ates asymptotically. In our scenario 
presented in section IV E the anomalously extended FDT 
regime is attributed to the soft droplets. It satisfies FDT 
but is absent in the ideal equilibrium where there are no 
frozen-in extended defects: it is a dynamical object. 
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FIG. 23: ZFC linear susceptibilities vs 1/L e (r) at T/J = 1.2 
and T/J = 0.8. The curves with small symbols are xzfc(i~ + 
tw,tw) with U, = 0, 10, 10 2 , 10 3 , 10 4 , 10 5 from the top to the 
bottom. The curves with solid lines are 1 — C(r + t w ,tw) of 
tw = 0, 10, 10 2 , 10 3 , 10 4 , 10 5 from the top to the bottom. The 
data of Txeq(r) (filled triangle) and its linear fit are shown 
at the bottom. The linear fit to T^zfc (t, 0) is also shown at 
the top. 



G. Integral Violation of FDT 

Let us finally examine the integral violation of the 
FDT defined in (||). In Fig. ||, the integral viola- 
tion I(t + t w ,t w ) obtained using our data is plotted 
against L(r)/L(t w ). We see it decreases with increas- 
ing waiting time t w deep in the quasi-equilibrium regime 
x = L(r)/L(t w ) <C 1 and increases deep in the aging 
regime x = L(r)/L(t w ) > 1. 

The intriguing problem is the validity of the FDT in 
the cr ossover regime x — L(r)/£(t w ) ~ 1. In section 
IV E 3 we conjectured that it is satisfied in the crossover 



regime. In Fig. 24, one can see that the integral violation 
is indeed decreasing with increasing t w atL(r)/L(t w ) ~ 1 
which appears compatible with our expectation. 

a possible scaling vari- 



As noted in section 1VE3 



+ 




1 1.5 

x=L(t)/L(tw) 



+ 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




10 2 
10 3 

10 4 
10 5 



X 




0.5 1 1.5 

x=L(z)/L(tw) 



FIG 



24: Integral violation of FDT at T/J = 1.2 (upper fig- 
and T/J = 0.8 (lower figure) plotted against L(r)/L(t w ) 



able for the interior of the crossover regime would be 



a = r/t w . In the limit t w — > oo, possible limits classified 
by different a are all smashed into the crossover regime 
x = L(r)/L(t w ) = 1. In Fig. ^5|, we show the integral 
violation against a = r/t w . Indeed, it is decreasing func- 
tion for any value of a = r/t w including even the highest 
a within our data. These observations are compatible 
with our conjecture that the FDT is asymptotically valid 
in the crossover regime. 



VIII. DISCUSSION AND SUMMARY 

To summarize, we have presented the results of a de- 
tailed MC simulation of a 4d EA Ising spin-glass model in 
the present paper. We demonstrated that the results can 
be well understood within the extended droplet theory 
we proposed recentlyEn 

We demonstrated the dynamical susceptibility xd 
larger than the equilibrium susceptibility xea exists in 
agreement with the extended droplet theory. Within the 
latter scenario, the two different limits emerge as dif- 
ferent large time (size) limits of the two-time (length) 
quantities. This phenomenon is presumably intimately 
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FIG. 25: Integral violation of FDT at T/J = 1.2 (upper fig 
ure) and T/J = 0.8 (lower figure) plotted against r/t w . 



related with the well known experimental observation 
Xfc > Xea which has been known, for a long time since 
the very discovery of spin glasseso. The question was, 
if it is a short time transient phenomena or not. The 
dynamical MFT was the first which clearly appreciated 
that the difference xfc — Xea (called "anomaly"El) exists 
in a well-defined large time limit. Our scaling theory pro- 
vides systematic extrapolation scheme to obtain the two 
susceptibilities in a controlled way. It will be certainly 
interesting to perform such an analysis in experiments 
and numerical simulations. 

We found evidences of the abrupt change of the two- 
time quantities in the crossover regime L(t) ~ L(t w ) 
which was anticipated in the extended droplet theory. 
This is directly related to the existence of the two differ- 
ent susceptibilities mentioned above. In previous experi- 
mental studies, it was well known that the relaxation rate 
S(t, t w ) has a broad peak as a function of t at around 
r ~ t w . This phenomenon is expected to be intimately 
related with the abrupt changes in the crossover regime. 
We proposed a modified relaxation rate S mo d(x, t w ) with 
x = L(t) I L(t w ). As expected, it was found numerically 
to develop a peak at x ~ 1 which sharpens with increas- 



ing t w when plotted against x. Probably the latter will 
be more useful in future studies. 

A very important feature of our analysis is the two- 
stroke strategy. In numerical simulations one can ob- 
tain directly the data of dynamical length scale L{t) by 
which time- dependent quantities are immediately trans- 
lated into length- dependent quantities. In addition, the 
information of the stiffness exponent 9 and fractal dimen- 
sion df of droplets are known from independent studies 
of equilibrium properties. Thus we had virtually zero fit- 
ting parameter left for us in this stage to test various 
data collapse expected from the scaling ansatz. 

The analysis of the growth law itself can be done sep- 
arately. We found that the growth law Lit) measured 
at various temperatures below T c shows the anticipated 
crossover from short time (length) critical dynamics to 
asymptotic activated dynamics. In our scaling analy- 
sis, we used the critical exponents v, z and the critical 
temperature T c determined in previous studies on equi- 
librium critical properties so that we have only one free- 
parameter ip which was found to be 2.5 — 3 in the present 
model. All the temperature dependences are renormal- 
ized into the crossover length Lq(T) and the correspond- 
ing crossover time scale ro(T) cx Lq(T) z by which all data 
are collapsed onto a universal growth law function. The 
latter becomes the expected power law for the critical 
dynamics at short times (length) and a slower function 
at large times (length) suggesting activated dynamics. 

Similar analysis should be done on experimental 
studiesEj and numerical simulations of other systems, 
particularly for 3d models, as well. In experiments, how- 
ever, this two-stroke strategy is not possible and proba- 



t quantities such 
to work out the 



bly the best way is to start from simp 
as relaxation of the AC susceptibility! 
parameters of the growth law. The advantage of the ex- 
periment, on the other hand, is that one can explore the 
order of 100 lattice spacings while the present numerical 
simulations are limited to 1 — 10 lattice spacings. 

Concerning the effects of critical fluctuation, many re- 
cent studies have pointed out the importance to renor- 
malize its effect. In three dimensions, previotts-^esults 
on the growth law Jay numerical simulationsE3EI'E2l as 
well as experimentsEj are well fitted to a power law as 
Lit) ~ i 1 / 2 ^) where the exponent l/z(T) is proportional 
to T, (but this T-linearity of X/ziT) starts, to break in 
the temperature range around T - Q.75TJB). Probably 
it is a sort of interpolation formula between T c , where 
critical fluctuations are dominant, and T — * 0, where all 
time scales associated with thermally activated processes 
diverge. Indeed a recent experiment in 3d systemsEj sug- 
gests the logarithmic growth law works well if critical 
fluctuation is properly considered. A recent numerical 
study of the growth lawE2l has also found a signature of 
the crossover. 

It is useful to note that the empirical power law 
Lit) ~ t 1 / z ( T ) j combined with the scaling laws in terms 
of L{t) 1 explains many of the empirical formulas pro- 
posed in previous experiments and numerical studies in a 
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unified manner being consistent with our two-stroke ap- 
proach. For example a fitting formula C oq (r) = ^ea + 
1/t q used in previous numerical studies (For example 
Ref. [30|) can be understood as a variant of ( |65|) with 
a = 8/z(T). The weU-Joiown "sub-aging" scaling for 
the TRM susceptibilityOEa xtrm(t + t w ,t w ) ~ F(t/<&) 
with \x < 1 can be understood as a variant of (j87]) wit i 
/i = 1 — 6/X. The "ojf -scaling" of the AC susceptibility!] 
can also be understood similarly as akeady noted in a 
numerical studja and an experimentalij study. In addi- 
tion, some apparent sub-aging feature of AC and DC sus- 
ceptibilities can be removed by considerations of finitc- 
ness of cooling rates in real experimentaij'Ed'E^I. 

The fundamental assumption of our scenario is that 
the effective stiffness constant T e fi of the free-energy gap 
FjJ of droplets is a function of the ratio x = L/R of the 
two length scales, namely the length scale of the droplet 
itself L and of the extended defect R which surrounds 
the droplet. Most importantly we assumed the vanishing 
of the effective stiffness constant as x — > 1 . This allowed 
the emergence of the two different order parameters <7ea 
and go and the associated susceptibilities xea and xd • It 
is desirable to clarify the scaling of the stiffness constant 
explicitly in the whole range of < x < 1. 

We expect our conjecture is consistent with the results 
of recent studies of low energy excitations in spin-glass 
modelala. In spin glasses of finite sizes R, it is likely 
that the existence of boundaries will intrinsically, induc e 
certain defects as compared with infinite systemsE3c£lcZl. 
Then we expect droplet excitations as large as the system 
size itself L ~ R is anomalously soft. Such an. anomaly is 
indeed found in the series of recent studieap. Although 
there new exponent 6' = was conjecturedli3, we consider 
it is better to attribute this to the zero stiffness constant 
T c ff[l] = as in (|43|). The stiffness exponent 9 > 0, on 
the other hand, is associated with a defect in T as we 
adopt in ([i"l|). The anomalously low energy and large 
scale excitations at L ~ R explains the apparently non- 
trivial overlap distribution function P(q) fo und .in nu- 
merous numerical studies of finite size systemsc§c3 which 
appear very similar to the prediction of the equilibrium 
mean-field theorjo. Although the meaning of the appar- 
ently non-trivial (and probaMy non-self averaging) P(q) 
in equilibrium is not obviousES, we consider the contribu- 
tion of the anomalous excited states to the macroscopic 
magnetic susceptibility in the present dynamical situa- 
tion, which is the realistic situatio n, is very important. 



As we discussed in section [VE 7, our scenario implies 



the average overlap q — dqqP(q) measured in equi- 
librium of finite size systems R (with built-in defects) 
is equivalent to the dynamical order parameter <7d- We 
found the latter is related to the Field Cooled suscepti- 
bility xfc as xfc = Xd = (1 - qD)/k B T (See @ and 
( |S6| ) ) . It would remind one of a folklore found in some lit- 
eratures that the difference of xfc and xzfc is somewhat 
related to difference of q and <7ea which is a consequence 
of the non-trivial P{q) found in Parisi's Miclica symmetry 
broken solution of the mean- field modeE3. However, we 
still have to recall the intrinsically dynamical nature of 
the situation: the domain walls, which allow the anoma- 
lous softening of droplets, are dynamical objects and they 
should be absent in ideal equilibrium where P(q) will 
have only one delta peak at q — qEA as predicted by the 
original droplet theoryU. 

Concerning the dynamical MFT, a serious problem in 
practice is that correction terms to the asymptotic limit is 
not known. They should obviously exist in the C— x rela- 
tion since the curves in Fig. ^systematically moves with 
increasing t w . Such a feature exist in the data of a re- 
cent experiment of simultaneous magnetic noisew 
measurements and previous numerical studiesH 

Our numerical results indeed suggests the separation of 
the breaking of TTI and FDT (see Fig. ||) .^nch^iea.- 
ture has not been realized in previous studicsHOiSEJEil. 
Within the scaling theory, the correction terms to the ex- 
pected asymptotic limit (Fig. || and Fig. |3|) themselves 
are predicted to have salient universal scaling properties 
which are amenable to be examined in practice. Our nu- 
merical results were well explained by the latter sca ling 
ansatz including the decay of the TRM susceptibility ( |87| ) 
which itself was predicted by the original droplet theory 
more than a decade ago. 
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